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SECTION  1.  INTRODUCTION 


For  a  structure  to  successfully  satisfy  the  needs  for  which  it  is  designed,  the  structural 
analyst/designer  is  often  required  to  predict  its  response  under  a  variety  of  dynamic  loading 
conditions.  Such  a  capability  to  predict  the  response  of  a  large  complex  structure  hinges  about  our 
ability  to  create  suitable  models,  and  to  be  able  to  improve  and  validate  them  through 
experimentatiqn  and  testing.  This  process  of  the  creation  of  an  analytical  model  for  a  structural 
system  along  with  procedures  to  experimentally  test,  upgrade  and  validate  in  some  sense  the  model 
is  what  shall  be  referred  to  in  this  report  by  the  term  structural  identification.  Thus  the 
identification  process  may  be  seen  as  a  necessary  prelude  to  the  active  control  of  structures  -- 
control  required,  for  instance,  to  suppress  vibrations,  achieve  fine  pointing  in  astronomical  space 
telescopes,  maintain  precise  shape,  limit  structural  stresses  caused  by  certain  types  of  dynamic 
loads,  or  satisfy  certain  tolerance  requirements  for  the  adequate  operation  of  certain  subcomponents 
and  precision  instrumentation. 

As  such  the  reason  why  the  concept  of  structural  identification  is  necessary  and  is  gaining 
increased  importance  as  we  continue  to  build  large,  complex,  one-of-a-kind  structures  is  our 
fundamental  inability  to  analytically  model,  with  desired  precision  and  confidence,  the  structural 
dynamics  of  large,  especially  flexible,  structures.  The  phases  of  model  building,  experimentation 
design,  algorithm  development  and  validation  procedures  are  actually  interactive,  each  dependent 
on  the  other,  and  all  in  turn  dependent  on  the  final  needs  which  the  desired  structure  is  meant  to 
satisfy.  This  interactive  nature  of  the  structural  identification  process  therefore  requires 
experimental  design  to  be  fully  cognizant  of:  (a)  the  analytical  model  to  be  used  (e.g.,  linear  or 
nonlinear  models),  (b)  the  nature  of  the  estimation  algorithm  to  be  used,  (c)  the  validation  criteria. 
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and  (d)  the  overall  goals  and  needs  that  the  structure  is  designed  to  meet. 

This  report  deals  with  the  development  of  a  methodology  for  the  optimal  placement  of 
sensors  in  structural  systems  for  system  parameter  identification^ PID)  and  some  associated 
aspects  of  this  general  problem  area.  It  may  therefore  be  thought  of  as  being  primarily  related  to  the 
experimental  design  aspects  of  the  SPID  process.  Yet,  as  we  shall  see,  the  interaction  with  the  other 
three  phases  of  the  process,  namely,  algorithm  development,  validation  and  model  building  is 
intense  and  cannot  be  neglected. 

1.1.  Motivation. 

Millions  of  dollars  are  being  spent  in  various  industries  related  to  the  fields  of  aerospace, 
mechanical,  civil  and  marine  engineering  for  the  experimental  testing  and  validation  of  structural 
systems  in  order  to  ascertain  that  they  behave  under  dynamic  loads  in  a  manner  that  they  are 
supposedly  designed  to.  Such  experimental  procedures,  in  essence,  are  rather  simple.  They 
constitute  the  application  of  known  dynamic  forces  to  the  structural  system  and  the  measurement  of 
its  response  at  one  or  more  locations.  Yet  the  procedure  is  extremely  expensive  especially  when 
carried  out  on  full-scale  systems,  such  as,  spacecrafts  and  building  structures.  To  gain  the  type  of 
information  that  would  provide  the  best  possible  model  for  a  structural  system  (and  thereby  provide 
the  capability  to  predict  its  response  to  a  variety  of  dynamic  forces  that  the  system  may  be  subjected 
to  during  its  useful  life)  one  therefore  needs  to  design  these  experiments  appropriately  so  that  a 
maximum  amount  of  information  can  be  obtained  about  the  relatively  unknown  aspects  of  the 
model  through  the  recorded  force  and  response  time  histories  of  motion.  Though  this  aspect  of 
experimental  design  has  been  long  recognized  by  experimentalists  and  structural  analysts ,  little,  if 


any,  work  has  so  far  been  done  in  this  area.  In  fact  a  literature  search  shows  but  a  few  journal 
references[  1  ]  on  developing  a  methodology  for  optimally  locating  sensors  in  a  dynamic  structural 
system  so  that  records  obtained  from  those  locations  would  be  maximally  beneficial  in  identifying 
the  structural  system. 

For  large  flexible  structures,  which  are  extended  in  space,  it  may  not  be  altogether  easy  to 
ascertain  where  to  locate  sensors  to  obtain  maximal  information  from  the  records  obtained.  The 
decision  on  where  to  locate  sensors  during  dynamic  testing  of  a  structure  is  generally  done  today 
by  experimentalists  and  technicians,  and  is  based  on  experience,  intuition  and  convenience  in 
acquiring  data.  It  is  unfortunately  not  based  on  any  cogent  theoretical  analysis  of  the  situation.  It  is 
such  an  analysis  that  this  report  attempts  to  include.  As  we  shall  show  in  this  report,  the  optimal 
sensor  locations  may  be  rather  nonintuitive  to  see  -  a  consequence  of  our  lack  of  experience  with 
inverse  problems.  Without  any  well  defined  methodology  for  locating  sensors,  so  that  data  from 
those  locations  will  be  maximally  useful  in  system  identification,  often,  experiments  which  yield 
little  or  no  useful  information  about  those  aspects  of  the  structural  model  that  the  analyst  is  greatest 
in  doubt  of,  are  carried  out  This  in  turn  necessitates  more  experiments,  all  adding  up  to  increased 
costs  in  terms  of  financial  and  human  resources,  each  new  experiment  being  designed  on  the  basis 
of  experience  and  intuition  rather  than  on  some  discernible  rational  methodology.  Such  an  ad  hoc 
experimental  design  process  is  not  only  extremely  wasteful,  but  does  not  ensure  that  we  are  doing 
the  best  we  possibly  can.  No  confidence  is  provided  in  the  results  from  a  particular  test  in  terms  of 
its  information  content  vis-a-vis  an  unknown  model  parameter  which  the  test  aims  at  assessing. 
Thus  not  only  would  such  testing  be  uneconomical  in  terms  of  time  and  money  but  it  would  lead  to 
erroneous  results,  and  worse  still,  a  lack  of  knowledge  regarding  how  badly  off  (or  erroneous) 


This  study  has  been  motivated  by  the  following  types  of  questions,  none  of  which  so  far 
have  been  adequately  answered  (several  have  not  even  been  addressed)  in  the  literature,  on  a 
rational  basis. 

•  Given  M  sensors,  where  should  they  be  placed  in  a  spatially  extended  dynamic 
structural  system  so  that  the  unknown  system  parameters  can  be  "best"  identified  ? 

•  Given  that  M  sensors  are  in  place  in  a  structure,  where  should  the  next  N  be  placed  so 
that  the  M+N  sensors  will  yield  maximal  information  about  the  unknown  system  parameters  ? 

•  What  is  the  minimum  number  of  sensors  required  in  a  given  structural  configuration  to 
be  able  to  uniquely  identify  a  certain  set  of  unknown  parameters  ? 

•  How  does  one  compress  data  in  the  recorded  time  series  during  an  experiment  without 
losing  information  about  the  parameters  which  need  to  be  identified  in  a  dynamic  system  ? 

•  Are  there  any  trade  offs  between  identification  and  control  in  so  far  as  the  placement  of 
sensors  is  concerned  ?  Are  there  sensor  locations  that  would  be  maximally  beneficial  for  control 
while  being  minimally  beneficial  for  identification  ? 

1.2.  Brief  Outline  of  the  Report. 

While  it  has  been  generally  recognized  that  the  quality  of  parametric  identification  of 
structural  systems  depends  on  the  location  at  which  sensors  are  placed  and  data  gathered,  very  little 
by  way  of  a  rational  methodology  has  so  far  been  developed  for  optimally  locating  sensors.  This 
report  presents  some  new  results  in  this  direction  and  points  out  the  need  for  more  work  —  work 
that  is  sorely  needed  if  the  costly  phase  of  experimentation  is  to  yield  quality  data  which  can  be 


used  to  best  identify  structural  systems.  While  we  do  not  pretend  that  we  have  entirely  solved  the 
five  questions  posed  in  the  previous  section,  several  new  results  pertaining  to  them  have  been 
obtained.  In  what  follows  we  provide  a  bird's  eye  view  of  the  topics  considered  in  each  section  of 
this  report,  sacrificing  at  times  exactitude  and  detail  (which  may  be  abundantly  found  in  the  material 
contained  in  the  individual  sections  themselves),  to  provide  a  general  sense  of  flow  and  direction  to 
the  entire  report. 

Section  II  begins  with  a  study  of  uniqueness  in  the  identification  of  a  specific  structural 
system  which  is  commonly  met  with  in  aircraft  and  space-craft  structures  and  more  generally  in  the 
fields  of  aerospace,  mechanical,  civil  and  marine  engineering  --  the  bending  cantilever  beam.  Thus 
we  look  here  at  situations  akin  to  those  posed  by  the  third  question  in  Section  1.1.  The  beam  is 
modelled  in  the  standard  manner,  using  finite  elements  and  its  mass  distribution  is  assumed  to  be 
known.  It  is  shown  that  for  such  an  undamped  bending  beam  model,  the  stiffness  distribution 
throughout  the  beam  can  be  uniquely  determined  through  the  placement  of  only  two  sensors  when 
the  beam  is  subjected  to  a  base  motion.  In  fact  if  the  displacement  and  rotational  time  histories  of 
any  two  adjacent  nodes  is  available,  unique  identification  of  the  entire  stiffness  distribution  of  the 
beam  to  one  side  of  the  fixed  support  is  possible.  While  this  may  be  thought  of  as  an  extension  of 
work  done  in  the  past  on  shear  beams[2],  the  results  obtained  here  will  find  wider  applicability 
because  of  the  ubiquitous  nature  of  bending  beams,  especially  in  space  structures.  While  these 
results  are  predicated  on  the  basis  that  there  is  no  measurement  noise,  they  indicate  an  important 
fact,  namely  that  sufficient  information  to  uniquely  identify  the  stiffness  distribution  is  available 
through  the  use  of  just  two  sensors,  appropriately  placed.  Also,  less  than  two  would  lead  to 
non-unique  identification,  as  would  the  use  of  two  or  more  sensors  if  placed  improperly  along  the 
cantilever.  As  most  identification  algorithms  are  iterative,  convergence  of  such  schemes  to  a 


particular  parameter  value  may  not  guarantee  that  we  have  found  the  correct  system  parameter 
unless  we  start  off  with  our  initial  guess  sufficiently  close  to  the  unknown  value.  The  objective 
function  surface  has  several  local  minima  in  which  such  an  iterative  scheme  could  get  trapped.  The 
significance  of  these  uniqueness  results  lies  in  the  fact  that  they  show  that  one  would  converge  to 
the  correct  parameter  values  starting  from  anv  initial  guess  should  we  appropriately  locate  our 
sensors;  furthermore,  only  two  sensors  are  required.  The  novelty  in  this  result  lies  in  that  even  the 
forward  solution  of  this  problem  cannot  be  obtained  in  closed  form. 

Section  HI  deals  with  the  development  of  a  methodology  for  optimal  placement  of  sensors 
when  one  is  interested  in  obtaining  the  "best"  estimate  of  the  unknown  parameters  starting  from  a 
suitable  close-by  initial  guess  and  using  noisy  measurement  data.  Thus  this  section  deals  with  the 
local  optimization  problem  where  one  is  attempting  to  reach  the  minimum  of  an  objective  function 
having  started  off  from  a  close  enough  initial  guess.  Loosely  speaking,  we  attempt  to  answer  the 
following  question:  a)  Given  that  we  are  attempting  to  fit  noisy,  measured  data  to  our  model  and 
thereby  determine  the  requisite  unknown  parameters  of  the  model,  and  b)  that  we  are  hunting  for 
the  best  parameters  that  fit  the  data  within  a  small  neighborhood  of  our  a  priori  knowledge  of  the 
parameters,  c)  where  should  we  place  sensors  so  that  a  minimum  covariance  estimate  of  the 
parameters  can  be  obtained  ?  Clearly,  these  locations  would  be  partially  dictated  by  the  estimator 
and  the  estimation  algorithm  used  and  therefore  the  sensor  locations  would  be  estimator  dependent. 
To  get  around  this  major  difficulty,  we  decouple  the  estimation  problem  from  the  optimization 
problem  by  invoking  the  concept  of  an  efficient  estimator,  thereby  making  the  methodology 
independent  of  the  estimator  used.  Thus  the  methodology  which  is  developed  answers  questions 
akin  to  the  fist  two  posed  in  Section  1.1.  In  fact,  as  we  shall  see,  it  provides  a  very  comprehensive 
picture  on  optimal  sensor  locations  for  parametric  identification  in  a  noisy  measurement 


environment. 


Section  IV  deals  with  an  important  aspect  of  experimental  design.  It  attempts  to  throw  light 
on  the  ways  in  which  one  can  minimize  measurement  data  handling  and  processing  in  SPID 
without  relinquishing  information  contained  in  noisy  data  about  the  unknown  parameters  which  one 
needs  to  identify  from  the  measurements.  The  study  is  of  a  preliminary  nature  and  deals  with  a 
simple  single  degree-of-freedom  system.  It  has  two  important  outcomes.  Firstly,  it  debunks  the 
myth  that  "the  larger  the  set  of  data  you  collect  (e.g.  frequency  data )  on  the  response  of  a  dynamic 
system  the  better  off  you  are  in  identifying  its  unknown  parameters."  In  fact,  we  analytically  show 
the  existence  of  frequency  data  that  would  have  no  bearing  on  reducing  the  uncertainty  of  the 
parameter  estimates  and  therefore  while  being  burdensome  in  handling  and  processing  would  cause 
no  improvement  in  the  results  of  SPID.  Such  studies  are  rudimentary  at  this  point,  but  will  indeed 
gain  importance  in  areas  such  as  in-flight  SPED  where  it  may  be  important  to  utilize  data  efficiently 
or  in  circumstances  where  data  needs  to  be  beamed  down  for  terrestrial  processing. 

Section  V  deals  with  the  tradeoffs  between  identification  and  control.  This  is  an  important 
aspect  of  SPID  as  it  is  subsumed  throughout  this  report  that  one  of  the  major  reasons  for  improving 
our  SPID  capability  is  so  that  we  can  control  large,  flexible  structures  to  perform  in  ways  that  are 
desirable  and  perhaps  meet  certain  specific  needs.  While  several  investigators  have  alluded  to  the 
duality  between  identification  and  control  on  a  qualitative  basis,  few,  if  any,  have  provided  any 
quantitative  analysis  of  the  situation.  In  this  section  we  do  this  through  the  intermediary  concept  of 
an  optimal  input  We  investigate  the  following  question:  Given  that  we  have  only  a  limited  amount 
of  energy  to  create  a  forcing  function  to  input  into  a  dynamic  system,  and  we  create  the  optimal 
forcing  function(s)  to  identify/control  the  system,  what  are  the  tradeoffs  between  our  ability  to 


identify  the  system  as  opposed  to  our  ability  to  control  the  system.  The  method  relies  on  the 
construction  of  an  objective  function  which  changes  continuously  from  one  that  emphasizes  control 
to  one  that  emphasizes  identification.  The  formulation  for  a  general  linear  dynamic  structural  system 
results  in  the  setting  up  of  a  nonlinear  two  point  boundary  value  problem.  While  the  method 
developed  here  is  applicable  to  structures  modelled  as  general  linear  dynamic  systems  we  have 
provided  in  this  report  numerical  results  for  only  a  single  degree-of-freedom  system. 

In  all  the  abovementioned  four  areas,  a  lot  remains  to  be  done.  In  fact  it  may  be  fair  to  say 
that  this  effort  barely  scratches  the  surface  of  the  deep  seated  problems  that  need  to  be  faced, 
identified  and  solved  if  SPID  of  large,  flexible  structures  can  be  carried  out  in  a  fruitful  and  reliable 
manner.  Yet,  each  section  of  this  report  contains  several  new  and  useful  results. 

In  years  to  come  the  role  of  SPID  will  indeed  increase;  for  the  determination  of  models  of 
complex,  flexible  structures  will  become  an  important  aspect  of  diverse  application  areas  like 
control  of  structures,  damage  assessment  of  structures,  early  fault  detection,  and  the  monitoring  of 
changes  in  structural  systems  either  as  a  consequence  of  environmental  conditions  or  as  a 
consequence  of  the  addition  or  removal  of  subassemblages.  An  excellent  paradigm  to  keep  in  mind 
would  be  spacecraft  structures,  which  need  to  be  controlled,  whose  configurations  may  be  greatly 
changed  through  docking  operations  with  other  space  structures,  whose  materials  may  degrade  in 
time  due  to  say  photon  radiation,  and  for  which  early  fault  detection  may  be  highly  beneficial. 

While  several  of  the  chapters  deal  with  general  linear  dynamic  systems,  and  develop 
generally  applicable  methodologies  to  answer  questions  such  as  those  posed  in  Section  1.1  for  such 
systems,  constraints  of  time  and  resources  have  caused  us  to  limit,  for  the  purposes  of  this  study, 


our  numerical  studies  to  only  the  simplest  numerical  examples.  Yet  these  simple  problems  for 
which  numerical  results  have  been  presented  in  this  report  are  perhaps  quintessential  to 
understanding  the  physics  of  the  phenomena,  without  which  insight  into  the  behavior  of  complex 
systems  would  be  hard  to  acquire.  More  work  is  clearly  called  for.  Hopefully,  the  results  reported 
herein  will  form  a  basis  for  perhaps  a  more  protracted  and  deeper  follow-on  study  which  could 
include  many  of  the  major  issues  dealt  with  here.  Accordingly,  Chapter  VI  addresses  the  issues  of 
what  this  study  has  taught  us,  in  which  directions  further  work  and  analyses  are  required,  what 
new  areas  have  been  uncovered  and  what  an  integrated  research  plan  addressing  some  of  these 
issues  should  endeavor  to  pursue. 
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SECTION  2.  UNIQUENESS  IN  THE  IDENTIFICATION  OF  A 

CANTILEVER  BENDING  BEAM 


2.1.  Introduction 

The  identification  of  large  structural  systems  is  an  area  of  research  that  has,  over  the  last  few 
years,  continued  to  attract  a  significant  amount  of  interest  all  over  the  world.  Most  such  methods  of 
parametric  identification  rely  on  iterative  algorithms  which  endeavor  to  successively  reduce  the 
mismatch  between  the  measured  responses  at  one  or  more  points  in  the  physical  system,  and  those 
predicted  by  a  system  model  which  purports  to  adequately  represent  the  physical  structure.  For  a 
given  time  history  of  measured  inputs,  applied  at  suitable  locations  in  the  structure,  the  system 
model  parameters  are  adjusted  (usually  automatically)  so  that  the  time  history  of  the  mismatch  is 
reduced  as  far  as  is  computationally  possible.  While  such  a  time  history  matching  scheme  will 
often  converge  to  a  set  of  model  parameters,  there  is  no  guarantee,  in  general,  that  the  parameters 
so  obtained,  correctly  represent  the  actual  physical  system  under  consideration.  The  "zeroing-in" 
of  the  iterative  scheme  to  an  erroneous,  albeit  physically  reasonable,  set  of  parameters  has  great 
importance  from  a  practical  standpoint;  for,  this  could  lead  to  substantial  errors  in  the  calculation  of 
quantities  of  engineering  significance  like  bending  moments  and  shear  forces  within  the  structure 

[31- 


In  a  previous  study  [2]  it  was  shown  that  a  structure  modelled  as  a  discrete  cantilever  shear 
beam  subjected  to  excitation  at  its  fixed  end  (base)  exhibits  a  degree  of  nonuniqueness  in  its 
identified  stiffness  distribution  which  is  dependent  on  the  location  at  which  the  time  history  of  the 
response  is  gathered,  and  subsequently  used  for  history  matching.  For  instance,  matching  of  the 


displacement  response  at  the  mass  level  closest  to  fixed  end  (first  mass  level)  was  shown 
to  cause  the  iterative  schemes  to  zero-in  on  the  correct  stiffness  distribution  throughout  the 
structure.  On  the  other  hand,  history  matching  of  the  response  gathered  at  the  N*  mass  level  of 
an  N  mass  system  could  lead  to  N!  different  stiffness  distributions,  all  of  which  would  be 
indistinguishable  in  so  far  as  their  "base  input  -  N*  mass  response"  pairs  are  concerned.  This 
realization  has  indeed  had  a  significant  impact  on  the  experimental  testing  of  spatially  extended 
structural  systems;  for  the  placement  of  sensors  can  critically  affect  our  confidence  in  the  identified 
parameters  obtained  through  history  matching  techniques. 

Whereas  several  cantilevered  structural  systems  used  in  aerospace,  civil  and  mechanical 
engineering  can  be  modelled  as  shear  beams,  several  others  may  be  more  appropriately  modelled  by 
bending  beams.  Cantilevered  bending  beams  are  widely  used  in  buildings,  bridges,  peers, 
offshore  structures,  mechanical  components  and  subassemblages,  marine  vessels,  nuclear  reactors, 
aircraft  and  spacecraft  structures  and  in  geotechnical  engineering.  Hence  such  beams  whose 
cross-sections  and/or  material  properties,  in  general,  vary  along  their  lengths  constitute  structural 
systems  that  are  generic  in  character  and  therefore,  from  a  practical  standpoint,  they  are  extremely 
important  This  section  concerns  itself  with  this  category  of  structureal  systems. 
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We  consider  a  finite  element  model  (FEM)  of  a  cantilever  undamped  bending  beam,  and 
investigate  the  conditions  under  which  unique  identification  of  the  stiffness  of  each  finite  element 
throughout  the  length  of  the  beam  can  be  achieved.  The  mass  distribution  is  assumed  to  be  known 
a  priori  throughout  the  structure,  presumably  from  design  drawings.  In  this  paper  we  consider 
structural  systems  subjected  to  excitations  (displacement  and/or  rotational  inputs)  at  their  fixed 
ends.  It  is  shown  that  unique  identification  of  the  stiffness  parameters  involved  in  the  finite  element 
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model  can  be  done  by  using  data  obtained  at  the  node  level  closest  to  the  fixed  end.  Generalizations 
to  uniquely  identifying  some,  though  not  all,  of  the  stiffness  parameters  in  the  FEM  model  are  also 
provided.  The  wide  usage  of  cantilever  bending  beams  in  various  fields  of  engineering  makes 
these  results  obtained  here  applicable  to  a  wide  set  of  application  areas. 


2.2.  System  Model 


Consider  a  structural  system  modelled  as  a  linearly  elastic,  discrete,  2N  degree  of  freedom 
bending  beam  as  shown  in  Figure  1(a).  The  length  of  the  i^  finite  element  is  taken  to  be  L;  and  its 
flexural  characteristic  to  be  (EI)i  (where  E  refers  to  the  modulus  of  elasticity  and  I  to  the  moment  of 
inertia).  We  shall  assume  that  both  the  mass  and  the  rotary  inertia  at  each  node  are  known  apriori 
as  are  also  the  lengths  L;  throughout  the  model.  Using  the  sign  convention  shown  in  Figure  1(b)  the 
stiffness  matrix  for  the  i*  element  can  be  expressed  as  [4] 
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where  ai  =  2(EI)I/L13 .  The  dynamical  equations  of  motion  can  then  be  written  as 


A  ..  A  A 

M  v  +  K  v  =  f 


(2) 


where  M  =  diag(m1,  m^.-.m^  )  is  the  2N  x  2N  lumped  mass  matrix,  containing  the  mass  and 


Stiffness 


rotary  inertias,  fc  is  the  assembled  stiffness  matrix  and  f  is  the  vector  of  external  loads.  The 
2N-vector,  v,  is  defined  as 


v  =  [w1,e1,w2,e2,...wN,eN]T  0) 

where  WT,  0;  are  the  absolute  displacement  and  rotation  at  node  i.  The  elements  m^,  i  =  1,  2,...N 
correspond  to  the  rotary  moments  of  inertia.  Since  the  mass  matrix  is  positive  definite,  equation  (1) 
can  be  renormalized  to 
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with  mQ  =  m  j  =1. 

The  identification  procedure  requires  the  determination  of  ai(  i  =  1,  2,...,N  given  the  knowledge  of 
the  matrix  M,  0q(O  and  UQ(t).  We  shall  assume  that  the  system  start  from  rest  so  that  v(0)  =  v(0) 

=  0.  The  masses  nr  ,  i  =  1,2,...,2N  and  the  elements  Li?  i  =  1,  2,...,  N  will  be  assumed  to  be 
known  and  positive.  Physical  interpretation  of  the  results  will  be  provided  as  we  go  along. 

2.2.1.  Some  Useful  Results 

In  this  section  we  prove  some  results  that  will  be  used  subsequently.  The  first  two  lemmas 
follow  from  the  positive  definite  nature  of  the  matrix  Ks.  We  provide  them  here  for  the  sake  of 
completeness  and  also  for  setting  out  notation  which  we  will  make  use  of  later  on. 

Lemma  1:  Given  the  matrix  L  where 


in  which  D  is  invertable, 


det(L)  =det(D)det(A-BD-1C) 


(9) 


Proof:  By  block  Gaussian  elmination 
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so  that  dct(L)  =  det(W)det(Z).  Noting  (10),  Z  =  D  and  W  =  A-BD  lC  and  the  result  follows. 


Lemma  2:  The  matrix  K$  has  no  zero  eigenvalue.  In  fact 


det(K  )  =  3N,+1  n  |>o 

k-s  H^k-lH^k  j 


(ID 


Proof:  We  can  express  as 
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Using  Lemma  l.detK  =  det(D)  det(K  )= -  det(K  ,) 

1  “aN^OM 


where  K  $  x  is  the  (2N-2(s-l)-2)  x  (2N-2(s-l)-2)  matrix  obtained  by  deleting  the  last  two  rows  and 
columns  of  K,  corresponding  to  a  system  from  mass  j  to  m^j.  Using  (14)  successively,  the 


result  follows. 


Let  P(  ■  denote  the  matrix  formed  by  deleting  the  i^row  and  j*  column  of  [KS-AJ],  Let  q;  ■ 
=  det[Pu  ].  Let  P0  0  =  (K$->J)  and  R(A.)  =  det(PQ  q).  The  we  have  the  following  result. 


Lemma  3: 
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Lt  (qj/^)  =  0 
X— »*• 


(16) 


Proof:  From  the  definition  of  q^(X)  and  R(A.)  we  have, 


R(A.)  =  (-X)  2N’2s+2  +  T(-X)2N-2s+1  +  ... ,  (17) 

qit  i  (X)  =  (-X)2"'25^  +  (T-cu)  (-XP2*  +  ....  (18) 

qit j(X)  =  (-l)H-i  c.  i(X)2N'2s  +  a(-X)2N*2s'1  + ....  i  *  j  (19) 

where  T  =  trace(Ks),  [KJ  =  [c;  •],  and  a  is  a  constant  which  can  be  obtained.  Taking  the  limits 
the  result  follows. 

2.3.  Identification  of  stiffness  distribution  of  finite  element  model 

Using  relation  (5)  and  taking  Laplace  Transforms,  and  replacing  the  transform  variables  by 

iV5C  we  get 


(KS-\I)U(X)  =  G(X) 


(20) 


where  U(X)  and  G(X)  represent  transformed  quantities  and  I  is  the  identity  matrix.  Solving  now 
for  US(X)  and  0S  (X.)  we  get 


U(X)  + 


e(X)=  +  qt^)q2,2(^)] 

*  R(X) 


Denoting 


.  1  Ps(^> 

A - - - 


AS-1(X)-  as  U  t(X) 


,  s  =  1,2,.. ,,N 


1  %(*-) 

p  /i[  \  t  s  “  1,2,«.mN 

b*-iW-  a,  us  j(X) 


and  using  (7)  we  get 


6  3L,  et-1(X) 

<m2s..m2I.3)‘/2  (m2,.im2!-2)1K  U>1(X) 


B, .,W- 


3L,  ,  L,  9,-l») 


,  s  =  1,...,N  .  (26) 


When  the  Lt  Asl(X)  exists,  it  shall  be  denoted  by  A^r  Similarly,  Bs  l  =  Lt  Bs 

*■-*»  X-x» 


^  ■  ■•  ■vvva  -V .v  a  aw 


v.v.v.AV.wAv.v.v.y 
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when  the  limit  exists.  Clearly  the  above  limits  exist  if  and  only  if 


Lt 


6,  ■!») 


(27) 


Lemma  4:  If  As 4  and  B$  4  exist,  then  they  cannot  both  be  zero. 

Proof:  Using  the  definitions  of  As  l  and  Bs  l  the  result  follows. 

Lemma  5:  Assume  that  As  4  and  Bs  x  exist,  and  m^ ,  m^  j,  Ls+1  >  0.  Then  at  least  one  of  the 
following  two  conditions  must  be  satisfied: 


B.  .  2  I  ^ 


,1/2 


(a)  A  *  0  and  -5'1  *  — 
A  L 


A,., 


(28) 


A  .  2 

s-1  ^ 


(b)  B  ,  *0  and  *  -L  . 
sl  B  3  »+1 

S-1 


H^.l 


1/2 


i,  ”2, 


(29) 


Proof:  We  show  that  when  condition  (a)  is  not  satisfied  condition  (b)  is  always  satisfied. 

(i)  Let  As  l=0;  then  condition  (a)  is  not  satisfied.  But  by  Lemma  4,  Bs_,*0also  As  l  / 
Bs4  =0*2/3  Ls+i(m2s.i/m2s)1/2  and  condition  (b)  is  satisfied. 


(ii)  Let  As_,  *  0  but  Bsl  /  As  l  =  2/L^j  (m2s/m2j.1)1/2.  Again  condition  (a)  is  not 


rr- 


^  v"  *.  *V 


satisfied.  But  this  implies  Bs  l  *  0  and  Bs_j  /  As  j  *  3/2Ls+1(m2s/m2s.1)I/2,  so  that  condition  (b) 
is  satisfied.  Similarly,  we  can  show  that  if  condition  (b)  is  not  satisfied,  condition  (a)  will  always 
be  satisfied. 


Lemma  6:  Given  nij,  j  =  2s-3,...2s,  and  L- ,  j  =  s,  s-*-l  and  Uj(X),  0j(X),  j  =  s-1,  s;  VA.,  the 
parameters  as  and  as+1  can  be  uniquely  identified  under  the  following  conditions 


(1)  As  j  and  Bs  j  exist 


(2)  Either 


(a,  \,,Oa„d|^*J-(^) 


or  (b)  B 


A,  i  2  f 

.  *  0  and  ■= —  *  —  L  ,  ~ 

-i  3  J 


Proof:  We  first  observe  (Lemma  5)  that  when  condition  (1)  above  is  satisfied,  at  least  one  of  the 
two  conditions  2(a),  2(b)  must  always  be  satisfied.  Using  (21)  to  (24),  we  have 


U(X) 


Usl(X)  R(X) 


LV.v .. 


V<V/V  a 

- - - (As  l^U  +  BS  Iq2  2)  • 

Usl(X)  R(X)  2  5  2 


We  divide  the  proof  into  two  steps. 
Step  1.  Determination  of  as 
If  As.j  *  0,  then  by  (33) 


XU(X) 


Ki  LAS1(X)U*-1(?L) 

R(X) 


Bs-1  ^2,1 
As-1  *1.1 


Noting  Lemma  3  we  have 


1  XU(X) 

a.  =  -  «—  Lt  — - - 

As-i 


Similarly,  if  Bs  t  *  0,  using  (34)  and  the  results  of  Lemma  3, 


1  X0s(X) 

a  *  = —  Lt  — - - 

B,-i  u..,<X) 


Note  by  Lemma  4  that  both  A  .  and  B  ,  cannot  simultaneously  be  zero. 


Step  2.  Determination  of  a 


_  *<  •  .  ■ 


s-t-1 

(i)  Condition  2(a)  satisfied: 

Expanding  [KS-XI]  in  terms  of  the  cofactors  of  the  first  column  we  have. 


Rft)  (6as  +  6as.l 

qu(X)  “2,-1 


(3as+iL,.f3aLs) 


1/2 

m^m 


1/2 
2s- 1 


6a 


s+1 


1/2 


**1.1 


^s+l^i+l  ^4.1 

(m2,+2m2,-l)1/2  ^ 


A  6a  +  6a 
=  (- 


s+1 


”2M 


-X)  +  D(X)  . 


Also  by  (21)  and  (23) 


RO) 

qua) 


U.,(X) 
5  U  (X) 


+  B  , 

s-1 


By  Lemma  3,  Lt  D(X)  =  0.  Also,  for  As  j  *  0,  using  (33)  and  Lemma  3, 


Lt  [aB  , 

1  s  $-1 


U 


q2  1  B  1  As  lUs  |(^)  1 

- ±i]=  Lt  [a  _Ll_L! - hL  j 

U  q,  ,  X-x»  A  .  Him  q,  , 


s-l 


7  V 


■>  s  v 


B  j  aSflLSfl  ~asLs  1 
~3  \(m2sm2s-l>1/2/  ' 


Equating  the  right  hand  sides  of  equations  (38)  and  (39)  and  taking  limits  as  we  get 


2  L  .  B  . 

a  [— _ *1L  i±  ] 

1  1/2  1/2  -r  J 

“to-l  “2.  A,-l 


s  ,)  ,  u  ,(X) 

Li -  Lt  [UaA  ,  — - ] 

3  *  -1  U  (X) 


2  LB. 

a.  [  — - —  +  »  -1  ]. 

1  1/2  1/2  —  J 

”2,-1  As-1 


The  result  therefore  follows. 

(ii)  Condition  2(b)  satisfied: 

The  determination  of  as+1  is  similar  to  the  above.  We  therefore  only  sketch  the  proof  briefly. 
Expand  R(X)  in  terms  of  the  cofactors  of  the  second  column  of  the  matrix  [KS-XI]  and  divide  by 
qi2-  Also  using  (22)  and  (24)  we  get 


q22(A) 


U,(X)  q,  (k) 


0S(X)  q22(X)  es(X) 


Equating  the  two  expressions  for  R(X)/q^  2(X)  and  taking  limits  as  we  get 


s+l  1/2  1/2  D 

^s-l  *1 


=  Lt  [X-aBsl  -i±l 

Ls+1  S  Sl  0 


aL 

S  S 

M 

\< 

1 _ 

3 

2L 

+  _ L 

Ls.l 

B,1 

1/2 

“2*1 

1/2 

. 

and  hence  the  result. 


Defining 


.■  A  1  Ps(X) 

A  .(X)  = - ^ - 


6  US.,(X) 


9-<*»  <*,< 


1/2  1/2 


B  ,0)  &  IJb^i  ■  -3L~  U-W  +  _ 


1/2  1/2 
2s  m2s-2 


for  s  =  1,  2....N  with  rr^  =  m  ,  =  1,  we  see  that 


exist  if  and  only  if 


(47) 


(48) 


(49) 


We  then  have  the  following  results. 

Lemma  7:  The  quantities  As  j,  Bsl,  As‘  - 1  and  Bs'  - 1  defined  in  equations  (25),  (26),  (44) 
and  (45)  are  such  that  at  least  one  of  the  following  two  conditions  is  always  satisfied. 

1.  As  l  and  B,.!  exist 

2.  As'  •  l  and  Bs'  - 1  exist  . 


es-l(X) 

Proof:  Condition  1  occurs  iff  Lt  -  £  v  ,  <  00  .  If  this  condition  is  not 


satisfied,  then  it  follows  that  ,Lt  - - —  <  y  .  <  «>,  and  therefore  condition  2 

0s  j(X) 


is  satisfied. 


Lemma  8:  Assume  that  As'  - 1  and  Bs'  - 1  exist  and  m^,  m^  j,  Ls+1  >  0. 
Then  at  least  one  of  the  following  two  conditions  must  always  be  satisfied 


t- 1  ^  2  ,  ”2*  ,1/2 


(a)  A  .  *  0  and  — *  r*-( — — ) 
-1  A'  m2s-. 


vb)  Bs, 


*  0  and  —  *  (^Li)1/2 


Proof:  The  proof  is  similar  to  that  of  Lemma  5. 


Lemma  9:  Given  mj,  j  =  2s-3,...2s;  Lj,  j  =  s,  s+1,  and  Uj(X),  0j(X),  j  =  s-1,  s;  for  all  X,  then 
parameters  as  and  a  S+1  can  be  uniquely  determined  under  the  following  conditions 


v  .■  / 


(1)  A‘  j  and  Bs'  - 1  exists,  and 


(2)  Either 


.  .  — '  _  ,  Bs-i  2  ,  An 

(a)  A  *  0  and  — 

A  .  S+1  “25-1 

5*1 


s*l  .2.1  /  “2S-1 


or  (b)  B  j  /Oand^i  *  f  L  (-^-)  . 

r  3  S+I  m.c 


Proof:  The  proof  is  along  the  same  lines  as  that  of  Lemma  4.  We  therefore  only  provide  the 


results. 


(i)  If  As'-i  *0 


a  =  -_  Lt  — i 


XU(X) 


=  [  2-..Hj±  hl±[  /^2lL  Lt  (X  +  a  — 


A  ,(X)1 


(ii)  If  Bs'-i  *  0 


a  =  — - —  Lt 

s  b\x- 


X8,(X) 


9S-,W 


,  and 


(57) 


a,+i“ 


[- 

1  / 

__L_  L  - 

3 

>1 

(/> 

1/2  S+1 
_m2s 

1/2 

“2.-1 

V 

t 

'  es 

aB  -ii - ] 

5  sl  6(\) 


a  L 

s  s 

X 

3 

+  < 

j 

LS.l 

™1/2 

HVl 

1/2 

m2s  _ 

j 

(58) 


We  note  that  at  any  node,  s-1,  of  the  finite  element  model,  (by  Lemma  7),  at  least  one  of  the 
conditions  (30)  or  (52)  must  always  be  satisfied;  thus  at  least  one  of  the  two  Lemmas  6  and  9  is 
always  applicable.  We  have  therefore  shown  that  if  the  displacement  and  rotational  time  histories 
of  motion  are  known  at  two  consecutive  nodes  of  a  FEM  of  a  bending  beam  along  with  the  mass 
and  rotary  inertia  properties  corresponding  to  those  two  nodes,  then  a  unique  identification  of  the 
material  property  of  the  element  in  between  the  nodes  and  also  that  of  the  element  above  the  upper 
node  can  be  obtained.  We  next  show  that  the  data  obtained  so  far,  uniquely  determines  the 
displacement  and  rotational  time  histories  at  the  next  node  point 

Lemma  10:  Given: 
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^  as’  as+l 

2.  Uj,  0J,  i  =  s-l,s 

3.  Lj,  i  =  s,  s+1,  and 

4.  tn^  i  =  2s- 1,  2s+2  , 


US+1(X)  and  0S+1(X)  can  be  uniquely  determined. 


Proof:  Considering  the  first  two  rows  of  equation  (20)  and  noting  the  structure  of  K  we  have 


u„,(W  = 


X  +  y  ^a, +1^1+1 

*  1/2 _ 1/2  *  1/2  1/2 


“2.  m2*+2 


“2,-1  “2,+2 


/Z,  • 


(59) 


and 


e»i<*)  = 


6a 


s+l 


*  1/2  1/2 


+  x. 


^as+l^s+l 


*  1/2  1/2 
“2s  m2s+l 


/z. 


(60) 


where. 


xs  =  Psa.)-(cu-x)us(X)-c12es(X) 


(61) 


Ys  —  qsft)  ■  ci^uj(X)  (c2(2  ^>)  0S(^) 


(62) 
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* 
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and 


z.  >° 


Hence  the  result. 


We  next  show  some  results  related  to  the  nature  of  the  quantities 

^s-l’  ®$-l’  ®s- 


Lemma  11.  If  As  t  and  Bs  l  exist  then 


(a)  if  As  *  0,  As  and  Bs  exist 


(b)  if  B^  0  ,  As'  and  Bs'  exist. 


Proof:  (a)  Using  (21)  and  (24), 


(65) 

U  A  ,q.  .  +  B  ,q,  , 

Dividing  the  numerator  and  denominator  by  VR(X),  and  taking  limits  X— and  using  Lemma  3, 
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we  get 


eft) 

Lt  s 

B  , 

s-l 

u,W 

A  . 

s-l 

Using  (25)  and  (26),  thus  As,  and  Bs  exist.  The  proof  of  part  (b)  is  similar. 
Lemma  12:  If  A  s  l  and  B's  l  exist  then 

(a)  if  A  s  l  *  0,  As  and  Bs  exist 

(b)  ifB's  l*0,As  and  Bs'  exist. 

Proof:  The  proof  is  similar  to  that  of  Lemma  1 1 . 

2.4.  Main  Result 

Theorem  1.  For  the  FEM  of  the  cantilever  bending  beam  considered,  given 

(1)  mj,  j  =  2s-3,...,  2N 

(2)  Lj,  j  =  s . N,  and 

(3)  Uj(t),ej(t),j*s-l,s, 


the  element  stiffness  distribution 


can  be  uniquely  determined. 


aj,j  =  s,  s+1 . N 


Proof:  The  proof  is  divided  into  three  parts  and  being  recursive  can  best  be  expressed  in 
algorithmic  form  as  follows. 


Do  steps  1  -*  3,  for  i  =  s,N-l 

Step  (1):  Find  A;  p  B; ,,  Ai  -  l  ,  Bi  -  t  using  expressions  (25),  (26),  (44) 
and  (45). 

Step  (2):  Use  either  Lemma  6  or  Lemma  9  depending  on  A;  j  and  Ai  -  l 
to  determine  a4  and  ai+1.  (By  Lemmas  7,  5  and  8,  one  of  these  two 
lemmas  must  always  be  applicable.) 

Step  (3):  Use  Lemma  10  to  find  Ui+1(X),  0i+1(X). 

We  have  thus  shown  that  the  stiffness  distribution  above  a  particular  node  of  the  FEM  can  be 
uniquely  determined  by  obtaining  the  displacement  and  rotational  time  histories  of  motion  at  that 
node  and  the  node  directly  above  it  Thus  the  distribution  can  be  uniquely  determined  by  locating 
sensors  at  only  two  points  in  the  system. 

Corollary:  The  stiffness  of  every  element  of  the  finite  element  model  of  the  cantilever  bending 
beam  can  be  determined  given 

(1)  mjf  j  =  1,  2,...,2N 

(2)  Lj,  j  =  1,...,  N,  and 

(3)  U0(t),  90(t),  Ujd),  0j(t). 

Proof:  Set  s  =  1  in  Theorem  1  and  the  result  follows. 
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if  the  base  displacement  and  base  rotational  time  histories  of  motion  are  known  as  well  as  the 
displacement  and  rotational  time  histories  of  motion  at  the  first  node  of  the  structure.  In  particular, 
unique  identification  of  the  entire  FEM  model  can  be  done  if  for  a  known  base  displacement  time 
history  of  motion  the  rotational  and  displacement  time  histories  of  motion  are  provided  at  the  first 
node. 

2.5.  Conclusions 

In  this  section  we  have  modelled  an  undamped  cantilever  bending  beam  by  the  finite  element 
method,  and  have  attempted  to  identify  the  stiffness  properties  of  the  elements.  The  stiffness 
matrix  of  each  element  is  taken  to  be  that  usually  used  in  day-to-day  structural  analysis  and  design. 
It  is  shown  that  knowledge  of  the  displacement  and  rotational  time  histones  of  motion  at  the  fixed 
end  and  at  the  node  closest  to  the  fixed  end  yield  unique  identification  of  the  properties  of  each  of 
the  finite  elements  of  the  model.  More  generally,  the  knowledge  of  the  rotational  and  displacement 
time  histories  of  motion  at  any  two  consecutive  nodes,  say  i  and  i+1  (see  Figure  1),  provides 
enough  information  to  uniquely  determine  the  stiffness  properties  of  all  the  elements  that  lie  to  one 
side  of  the  node  i,  away  from  the  fixed  end. 

It  should  be  noted  that  the  analysis  presented  here  deals  with  noise  free  data.  The  presence 
of  noise  in  the  measurements  would  in  general  lead  to  low  signal-to-noise  ratios  for  measurements 
at  the  lowest  node  levels  so  that  though  the  identification  problem  has  a  unique  solution,  in  actual 
practice,  the  variance  of  the  parameter  estimates  may  become  large.  Measurements  made  farther 


away  from  the  fixed  end  would  be  larger  in  amplitude  so  that  the  noise  related  estimation  errors 
would  be  smaller.  However,  one  faces  then  the  nonuniqueness  problem.  Furthermore,  the 
analysis  herein  assumes  that  a  complete  knowledge  of  the  input-output  time  histones  (or  their 
transforms)  is  available. 


While  no  attention  has  been  given  to  the  process  of  arriving  at  an  adequate  finite  element 
model  (i.  e.  the  proper  location  and  number  of  nodal  points)  the  prediction  of  the  response  of  the 
actual  physical  system  to  various  input  time  histories  using  the  identified  parameters  (obtained  by 
history  matching)  critically  depends  on  the  adequacy  of  the  model.  Structural  modelling  being  to  a 
good  extent  an  an  acquired  mainly  by  practice,  this  aspect  of  the  problem  has  not  been  considered. 


However,  it  should  be  pointed  out  that  the  uniqueness  results  arrived  at  here  are  ii 


i.  e.  no  matter  how  many  finite  elements  are  used  to  represent  the  system 


(i.  e.  no  matter  how  many  parameters  to  be  it 


(at  the  fixed 


end  and  at  the  node  nearest  the  fixed  end)  are  sufficient  to  tie  down  the  estimates  of  the  stiffness  of 


each  element  uniquely.  The  analysis  of  damped  structural  models  has  been  specifically  excluded 
for  the  sake  of  simplicity  and  work  along  these  lines  is  continuing. 


Cantilevered  bending  beams  form  one  of  the  commonest  categories  of  structural  elements  and 
are  used  widely  across  diverse  fields  of  engineering.  They  find  application  in  the  fields  of  civil, 
mechanical,  nuclear,  aerospace  and  marine  engineering  to  name  but  a  few  and  hence  these  results 
may  have  significance  to  a  wide  set  of  application  areas.  Large  boom  type  structures  in  the 
aerospace  industry,  tall  buildings  in  civil  engineering  and  the  identification  of  subassemblages  in 
the  nuclear  industry  are  but  a  few  specific  examples. 
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SECTION  3.  A  METHODOLOGY  FOR  OPTIMUM  SENSOR  LOCATIONS  FOR 


PARAMETRIC  IDENTIFICATION  OF  DYNAMIC  SYSTEMS 

3.1.  Introduction 

Reliable  predictions  of  structural  responses  are  closely  dependent  upon  the  validity  of  the 
models  chosen  to  represent  the  systems  involved.  When  parametric  models  are  used,  a  proper 
knowledge  of  the  various  parameter  values  becomes  crucial  in  establishing  the  usefulness  of  such 
models.  However,  to  actually  come  up  with  these  parameter  values,  one  often  needs  to  collect 
response  data  from  instruments  located  at  various  positions  within  the  structure.  The  usefulness  of 
such  data,  in  turn,  depends  primarily  upon  the  instrument  characteristics  and  on  the  chosen 
positions  where  the  instruments  are  located.  Consequently,  for  given  types  of  instruments,  which 
are  to  be  used,  one  often  wants  to  locate  them  such  that  data  collected  from  those  locations  yield  the 
"best"  estimates  of  the  modelled  structural  parameters. 

Although  various  methods  have  been  developed  to  identify  the  parameters  that  characterize 
flexible  structures  (e.  g.,  [5-10] ),  from  records  obtained  in  them  under  various  loading  conditions, 
few  investigators,  if  any,  have  looked  at  the  question  of  where  to  locate  sensors  in  a  large,  spatially 
extended  structure  to  acquire  data  for  "best"  parametric  identification  [1],  The  problem  of  optimally 
locating  sensors  in  a  dynamic  vibrating  system  mainly  arises  from  considerations  of:  (1) 
minimizing  the  cost  of  instrumentation,  data  processing  and  data  handling  through  the  use  of  a 
smaller  number  of  sensors,  data  channels,  etc.;  (2)  obtaining  better  (more  accurate)  estimates  of 
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model  parameters  from  noisy  measurement  data;  (3)  improving  structural  control  through  the  use  of 
superior  structural  models;  (4)  efficiently  determining  structural  properties  and  their  changes  with  a 
view  to  acquiring  improved  assessments  of  structural  integrity,  and  (5)  improving  the  early  fault 
detection  capability  for  large,  flexible  structural  systems. 

The  problem  addressed  in  this  section  can  succiently  be  stated  as  follows:  Given  m  sensors 
where  should  they  be  located  in  a  spatially  distributed  dynamic  system  so  that  records  obtained 
from  those  locations  yield  the  ’’best"  estimates  of  those  relatively  unknown  parameters  which  need 
to  be  identified? 

In  the  past,  the  optimum  sensor  location  problem  (OSLP)  was  solved  by  positioning  the 
given  sensors  in  the  system,  using  the  records  obtained  at  those  locations  with  a  specific  estimator, 
and  repeating  the  procedure  for  different  sensor  locations.  The  set  of  locations  which  yield  the 
"best"  parameter  estimates  would  then  be  selected  as  optimal.  The  estimates  obtained  of  course 
depend  upon  the  type  of  estimator  used.  Thus  the  optimal  locations  are  estimator  dependent,  and 
an  exhaustive  search  needs  to  be  performed  for  each  specific  estimator.  Such  a  procedure  besides 
being  highly  computationally  intensive  suffers  from  the  major  drawbacks  of  not  yielding  any 
physical  insight  into  why  certain  locations  are  preferable  to  others. 

Work  on  the  solution  of  the  OSLP  was  perhaps  first  done  by  Shah  and  Udwadia  [1].  In 
brief,  they  used  a  linear  relationship  between  small  perturbations  in  a  finite  dimensional 
representation  of  the  system  parameters  and  a  finite  sample  of  observations  of  the  system  time 
response.  The  error  in  the  parameter  estimates  was  minimized  yielding  the  optimal  locations.  In 
this  section,  we  developed  a  more  direct  approach  to  the  problem  which  is  both  computationally 


superior,  and  throws  considerable  light  on  the  rationale  behind  the  optimal  selection  process.  The 
methodology  is  applicable  to  all  spatially  extended  dynamic  systems.  In  this  report,  special  attention 
has  been  paid  to  large,  flexible  structural  systems  such  as  large  space  structures. 

We  uncouple  the  optimization  problem  from  the  identification  problem  using  the  concept  of 
an  efficient  estimator  (e.  g..  the  maximum  likelihood  estimator  as  the  time  history  of  data  becomes 
very  large)[  1 1],  For  such  an  estimator  the  covariance  of  the  parameter  estimates  is  a  minimum. 
Using  this  technique  and  motivated  by  heuristic  arguments,  a  rigorous  formulation  and  solution  of 
the  OSLP  is  presented. 

3.2.  Uniqueness  of  Identification  and  Local  Optimization 

The  iterative  nature  of  most  identification  schemes  require  us  to  differentiate  between  two 
distinct,  and  often  times  confusing,  criteria  for  the  sensor  location  problem. 

(1)  Uniqueness  of  Identification 

The  idea  here  is  to  locate  sensors  in  such  a  manner  that  no  matter  what  initial  guess  (of  the 
parameters  or  functionals)  one  starts  off  with  in  the  iterative  scheme,  the  identification  will 
converge  to  the  unique  "actual”  system.  Alternatively  put,  one  wants  to  locate  sensors  at  locations 
which  yield  information  that  can  unequivocally  tell  us  the  parameters  of  the  actual  system.  Using 
records  obtained  from  such  locations,  different  initial  guesses  would  yield  different  estimates  of 
the  parameters  (or  functionals)  from  those  of  the  actual  system.  Section  2  of  this  report  dealt  with 
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this  aspect  for  a  cantilever  bending  beam,  a  commonly  used  substructure  in  structural  design. 

Such  nonuniqueness  problems  if  they  in  fact  did  exist  could  indeed  lead  to  not  only  incorrect 
identification  of  structural  system  parameters,  irrespective  of  the  actual  estimator  used,  but  could 
also  lead  to  incorrect  estimates  of  quantities  of  engineering  significance  like  bending  moments  and 
shear  forces  within  the  structure.  Such  incorrect  estimates  would  directly  affect  our  ability  to  assess 
the  structural  integrity  of  the  system  and  our  ability  to  control  it  appropriately. 

(2)  Local  Optimization 

If  however,  one  has  a  fairly  good  idea  of  the  parameter  values,  the  initial  guesses  would  be 
fairly  close  to  the  actual  values.  Thus,  provided  the  measurements  are  not  too  noisy,  one  would 
not  be  likely  to  converge  to  a  solution  other  than  that  represented  by  the  actual  system.  Having 
restricted  the  search  space  in  this  manner,  so  to  speak,  through  the  use  of  a  good  initial  guess,  one 
needs  to  locate  sensors  in  such  a  manner  that  having  started  with  these  approximate  estimates,  the 
records  obtained  have  the  greatest  information  to  improve  the  estimates.  As  opposed  to  the 
"global"  convergence  (starting  from  "any"  initial  guess)  which  Section  2  addressed,  this  section 
looks  at  the  development  of  a  methodology  for  optimally  locating  sensors  in  the  context  of  "local" 
convergence  (starting  from  a  "close"  initial  guess).  In  Section  3.9.3  we  will  further  illustrate  some 
of  the  differences  between  these  two  concepts. 


3.3.  Model  Formulation 
3.3.1.  System  Model 

The  study  of  physical  systems  which  can  be  adequately  represented  by  linear  constant 
coefficient  differential  equations  will  be  considered  in  this  sequel.  Though  most  large  complex 
dynamic  systems  are  spatially  continuous  in  nature,  often,  suitable  discrete  models  can  be 
formulated  for  engineering  applications.  In  fact  finite  element  and  finite  difference  methods  have 
become  common  in  the  reduction  of  physically  continuous  systems  to  mathematically  discrete 
models. 

Though  the  development  of  the  optimum  sensor  location  (OSL)  criterion  will  be  shown  to  be 
unrelated  to  the  nature  (linear,  nonlinear,  time-variant,  time-invariant)  of  the  system  S  under 
consideration,  let  us  for  the  moment  consider  a  linear  dynamic  vibrating  system  so  that  we  have  a 
vehicle  for  developing  the  methodology. 

The  governing  differential  equation  of  morion  for  a  linear  dynamic  system  may  be  considered 

as  : 

MX  +  CX  +  KX  =  F(t) ;  X(0)  =  XQ  ,  X(0)  =  XQ  (69) 

where  X0  and  Xq  are  the  given  initial  conditions  for  the  system.  The  constant  coefficient  matrices 

M,  C  and  K  are  each  of  dimension  (N  x  N).  The  M  matrix  may  be  considered  as  the  lumped  or 
the  consistent  mass  matrix  (as  is  often  the  case  in  structural  analysis),  C  as  the  damping  matrix  and 
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K  as  the  stiffness  matrix.  X  is  an  N-vector  whose  components,  Xj,  may  be  considered  to  be  the 
displacement  response  (at  the  nodes  of  the  finite  element  or  the  finite  difference  mesh)  of  the  system 
shown  in  Figure  3  to  the  input  vector  {F(t)}N.  One  or  more  elements  of  the  coefficient  matrices, 
in  equation  (1),  may  constitute  the  unknown  parameters.  To  "best"  estimate  these  parameters,  one 
would  locate  sensors  in  the  system  in  such  a  way  that  the  measurements  obtained  thereat  are  most 
informative  about  the  estimated  parameters.  To  accomplish  this  task,  let  us  collect  all  the  possible 
unknown  parameters  in  a  vector  0  of  dimension  L.  Hence 

0  =  <  0M  I  0C  1  0K  >  T  (70) 

where  the  subvectors  0M,  0C,  and  0K  have  dimensions  a,  b,  and  c,  respectively.  The  superscript  T 
indicates  vector  transpose  throughout  this  dissertation.  In  general,  a,  b,  and  c  are  each  at  most 
equal  to  N2;  however,  in  many  dynamic  systems,  the  coefficient  matrices  in  equation  (69)  are 
symmetric.  In  such  a  case  a,  b,  and  c  are  each  at  most  equal  to  (N  +  l)N/2. 

3.3.2.  Measurement  Model 

The  response  of  the  dynamic  system  is  assumed  to  be  measured  using  an  m  (m  <  N) 
available  sensors.  Solution  of  the  OSLP  is  equivalent  to  the  selection  of  the  m  locations  out  of  N 
possible  such  locations  so  that  the  m  time  history  response  obtained  thereat  yield  the  maximum 
amount  of  information  about  the  system  parameters. 

To  formulate  the  measurement  model,  let  us  first  assume  there  are  exactly  N  sensors 


available,  so  that  each  component  of  X  is  measured  (i.  e.,  m  =  N).  These  measured  responses  can 
be  mathematically  represented  by  the  N-vector  Z  as  follows: 

Zj(t)  =  gj  [X(0,  t)  ]  +  Nj(t)],  j  =  1,2. . N  (71) 

where  Z ■  is  the  j-th  component  of  Z(t),  functional  gj  represents  the  "measurement  process",  and  the 
dependence  of  the  response  X  on  the  parameter  0  is  explicitly  noted.  We  shall  assume  that  g  is  a 
memoryless  transformation  of  the  system  output  which  yields  the  measurements.  The 
measurement  noise  Nj(t)  is  taken  as  nonstationary  Gaussian  white  noise  with  a  variance  of  \j/2(t). 
Therefore, 

EIN/tjJN/tj)]  =  V2(tj)  SK(i-j)  5D(t,  -  t*)  (72) 

where  5K  and  Sp  stand  for  the  Kroneker  and  the  dirac-delta  functions,  respectively.  Having 
measured  each  element  of  the  response  vector  X,  a  total  of  m  out  of  N  responses  need  to  be 
selected  so  that  they  contain  the  most  information  about  the  system  parameters  and  are  maximally 
sensitive  to  any  change  in  the  parameter  values.  This  "selection"  process  can  be  represented  by  an 
m-dimensional  vector  Y  such  that 

Y(t)  =  SZ(t)  (73) 

where  S  is  the  (m  x  N)  upper  triangular  selection  matrix  with  each  row  containing  null  elements 
except  for  one  which  is  unity.  The  m  different  components  of  Z  selected  to  be  measured  are  so 
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ordered  in  vector  Y,  that  if  the  element  in  the  i-th  row  and  k-th  column  of  S  is  unity,  the  (i+l)-th 
row  can  have  unity  in  its  i-th  column  only  if  i  >  k.  The  matrix  S  then  has  the  property  so  that 
matrix  P  =STS  is  an  (N  x  N)  diagonal  matrix  with  unity  in  its  i-th  row  if,  and  only  if,  Zi  is  selected 
to  be  measured.  The  elements  of  P  are  otherwise  zero.  Hence,  one  can  write 


Y(t)  =  Sg[X(9,  t)]  +  SN(t), 


Y(t)  =  H[X(0,  t)]  +  V(t) 


If  g;  is  linearly  related  to  the  response  x^,  in  general,  then 


H[X(0,  t)]  =  SRX 


where  R(t)  in  equation  (76)  can  be  thought  of  as  a  dynamic  gain  matrix.  In  the  case  that  g;  is 
related  to  the  response  only,  matrix  R  will  reduce  to  a  diagonal  matrix,  Diag(pj,p2,. . .,  pN). 

The  problem  of  locating  sensors  in  an  optimal  manner  then  reduces  to  determining  the 
selection  matrix  S  defined  above.  Alternately  put,  one  needs  to  determine  the  m  locations  along  the 
diagonal  of  the  matrix  P  that  should  be  unity.  These  locations  must  be  so  chosen  as  to  obtain  the 
"best"  parameter  estimates. 


3.4.  Efficient  Estimator  and  the  Cramer-Rao  Lower  Bound 
3.4.1.  Conditional  and  Unconditional  Estimation 

In  this  section  we  shall  consider  two  types  of  estimations.  (1)  The  conditional  estimation 
problem  results  when  0  is  an  unknown  constant,  and/or  when  the  a  priori  probability  density 
function  P(0),  is  not  known.  For  this  type  of  problem  the  expected  value  of  the  parameter  estimate 
equals  the  true  value  of  the  quantity  being  estimated  and  is  often  referred  to  as  conditionally 
unbiased  estimation  [12,13].  (2)  In  the  unconditional  estimation  problem  the  unknown  parameter 
0  is  characterized  by  a  known  prior  probability  density  function.  In  this  type  of  problem  the 
expected  value  of  the  estimate  equals  the  expected  value  of  the  quantity  being  estimated  and  is  called 
unconditionally  unbiased  estimation  [12,13]. 

Unconditional  estimation  is  more  general;  however,  in  many  situations  when  the  probability 
distributions  of  the  parameter  estimates  are  not  know,  one  can  not  utilize  them.  Suppose  it  is 
desired  to  find  an  estimate  of  the  parameter  vector  0;  if  the  joint  density  P0(0)  is  not  in  hand  or  0  is 
not  a  random  vector,  then  conditional  estimation  could  be  used. 

In  the  analysis  of  flexible  structures,  an  example  of  conditional  estimation  may  involve 
determination  of  the  structural  parameters  of  a  space  structure  from  noisy  observations.  In  such  a 
case,  the  unknown  parameters,  say  the  member  stiffnesses,  are  considered  as  unknown  constants 


which  remain  constant  even  after  many  observational  measurements  are  made.  An  example  of 
unconditional  estimation,  on  the  other  hand,  involves  some  knowledge  of  the  probability  density  of 
the  unknown  random  parameters. 


3.4.2.  Cramer-Rao  lower  bound  (CRLB)  and  Cramer-Rao  (CR)  inequalities 


The  lower  bound  of  the  covariance  of  the  estimation  error  was  first  proposed  by  Cramer  and 
is  called  the  Cramer-Rao  Lower  Bound  (CRLB)  [11,  14].  This  CRLB  is  related  to  the  minimum 
error  between  the  estimates  and  the  actual  value  of  the  parameters  to  be  identified  in  the  sample 
distributions  given.  The  CRLB  is  a  nontrivial  lower  bound  on  the  covariance  of  the  estimate. 
Unfortunately,  the  CRLB  depends  upon  the  bias  of  whatever  estimator  is  used.  However,  for 
unbiased  estimators  the  CRLB  is  indeed  independent  of  the  estimator  that  might  be  employed.  For 
any  chosen  arrangement  of  sensor  locations,  expressions  for  the  conditional  CRLB  (CCRLB)  and 
the  unconditional  CRLB  (UCRLB)  for  unbiased  estimators  are  derived  in  [15].  They  are: 


UCRLB  = 


P  (0)  •  (CCRLB)' 
0 


d0 


(78) 


Hence 
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E[  (0-0)  (0  -  0)T  1  0  ]  >  CCRLB 

(79) 

and 

E[(0  -  0)  (0  -  0)T,0]  £  UCRLB 

(80) 

The  right  hand  side  of  the  inequalities  (79)  and  (80)  can  be  obtained  from  equations  (77)  and  (78), 
their  left  hand  sides  depending  on  the  type  of  estimator  used.  Inequalities  (79)  and  (80)  are  known 
as  the  conditional  and  the  unconditional  estimation  covariance  inequalities. 

3.4.3.  Efficient  estimators 

If  an  unbiased  estimator  achieves  the  CRLB,  the  estimator  is  called  efficient.  An  efficient 
estimator  may  also  be  called  an  optimal  unbiased  estimator  since  it  achieves  a  minimum  "Expected 
Square  of  Estimation  Error"  (ESEE)  and  no  other  unbiased  estimator  could,  in  fact,  achieve  a  lower 
ESEE.  It  is  of  some  interest  to  also  note  (see  equations  79  and  80)  that  an  unbiased,  efficient 
estimator  is  also  a  minimum-variance  estimator.  It  is  however,  unfortunate,  that  there  are  no 
general  methods  for  constructing  efficient  estimators.  But  one  commonly  useful  class  of 
estimators  which  is  asymptotically  unbiased  and  efficient  is  the  one  which  comprises  the 
conditional  and  unconditional  maximum  likelihood  estimators  ([14],[16]).  Hence  for  efficient 
estimators  (minimum  covariance)  the  inequalities  [5]  become  equalities.  Therefore,  one  can  write 

E((0  -  0)  (0  -  0)T]  =  CRLB  (81) 

In  this  sequel  such  as  estimator  shall  be  assumed  to  exist. 
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The  expression  inside  the  brackets  in  equalities  (77)  and  (78)  are  known  as  the  Fisher 
Information  (F.  I.)  matrices  [15].  Therefore  maximization  of  this  Fisher  Information  matrix 
(maximizing  a  certain  norm  of  the  matrix;  such  as  the  trace  norm,  etc...)  would  yield  the  minimum 
possible  value  of  the  covariance  of  the  estimation  error  [17,18].  It  is  of  some  interest  to  note  that 
when  using  efficient  estimators,  the  covariance  of  the  estimation  error  is  known  with  no  regard  to 
the  actual  form  of  the  estimator. 


The  Fisher  Information  Matrix,  Q(T),  for  conditional  estimation  of  0  (using  equation  (77) ) 
can  now  be  written  as: 
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3.5.  Optimum  Sensor  Location  (OSL)  for  Conditional  estimation  ofVector-Valued 
Parameters 

3.5.1.  Some  Motivative  Thoughts  and  the  Fisher  Information  Matrix  (F.  I.) 


Consider  a  case  in  which  one  tries  to  estimate  one  parameter,  0j,  which  is  to  be  identified  in 
a  dynamic  system  model  with  only  one  sensor  provided.  One  would  want  to  ideally  choose  a 
location  i  (out  of  N  possible  such  locations)  such  that  the  measurement  y;  (t),  i  £  [l,  N],  t  e  (0,  T)  at 
location  "i"  yields  the  "best”  estimate  of  the  parameter  0r  Heuristically,  one  should  place  the 
sensor  at  such  a  location  that  the  time  history  of  the  measurements  obtained  at  that  location  is  most 


sensitive  to  any  changes  in  the  parameter  9j.  Hence,  in  equation  (75)  it  is  really  the  slope  of 
H[X(0j,  t))  with  respect  to  0j  that  needs  to  be  maximized.  However,  since  only  the  absolute 
magnitude  of  this  slope  is  of  interest,  it  is  logical  to  want  to  find  i  (or  equivalendy  determine  the 
selection  matrix  S  described  previously)  such  as  to  maximize  (3H/90j)2.  Since  this  quantity  is  a 
function  of  time  one  would  want  to  locate  a  sensor  which  maximizes  the  "average"  value  of 
(9H/90j)2  over  the  time  interval  (0,  T)  during  which  the  response  is  to  be  measured.  This  leads  to 
maximizing  the  following  integral: 

q,(T)  = 

When  there  is  more  than  one  parameter  to  be  estimated,  and  the  number  of  sensors  is  greater 
than  unity,  this  intuitive  approach  needs  to  be  extended  in  a  more  rigorous  manner.  In  such  cases 
recourse  to  mathematical  treatment  is  necessary,  and  we  shall  see  that  such  treatment  will  be  in 
agreement  with  our  heuristic  solution. 


3.5.2. 


Fisher  Information  Matrix 
vector-valued  problem 


in  terms  of  the  measurement  response  for 
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It  is  shown  that  for  any  conditionally  unbiased  estimator  of  0  the  covariance  estimation  error 
can  be  written  as 
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(84) 


El (0  -  0)(8-  $)T] 


The  integral  on  the  right  hand  side  of  the  inequality  (84)  is  the  well  known  Fisher  Information 
matrix,  as  defined  previously  in  equation  (82).  In  order  to  reduce  the  error  in  the  estimates,  one 
would  equivalently  want  to  maximize  a  suitable  norm  (e.  g.,  Trace,  etc.)  (16]  of  the  Fisher 
Information  matrix  Q(T).  Therefore,  introducing  equation  (76)  into  equation  (82),  (this  constitutes 
an  extension  of  the  equation  (83)  which  we  heuristically  derived  for  the  scalar  case,  to  the  vector 
situation),  one  obtains: 


xtrtprx 
— idt, 
V  (t) 


where  the  ij  clement  of  Xe  can  be  written  as: 


dx 

[X  L  =  — -  ,  i  e  [1,N],  j  e  (l,m) 

e  J  ae 

j 


(85) 


(86) 


where  X  =  {xt}  N  and  9  =  {0j}L.  We  note  that  the  Fisher  Matrix  is  symmetric  and  is  dependent  on 
the  length  of  the  record  available,  as  well  as  the  locations  of  the  sensors  as  determined  by  the  matrix 
P. 

If  the  m  locations  where  the  sensors  are  to  be  placed  are  denoted  by  sk,  k  =  1,  2,...,  m,  then 
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where  the  (N  x  N)  diagonal  matrix  1^  has  all  its  elements  equal  to  zero  except  the  element  of  the 
sk  row,  which  is  unity.  Noting  that  P  is  a  diagonal  matrix,  equation  (85)  can  be  simplified  to  yield 


m  J  XTrT  r  X  dt 

7  f  e  \  *k  e 


V  f  9 

';sps2 . sm;  S,0;I]  =  2*  J  - 


where  r^  is  the  sk  row  of  the  matrix  R.  Also  in  eq.  (87)  explicit  mention  is  made  of  the 
dependence  of  the  Fisher  Matrix  on  the  time  length  T  of  the  available  data,  the  system  S,  the 
parameter  vector  0,  and  the  time-variant  input  I.  If  the  matrix  R  is  diagonal,  with  diagonal  elements 
Pi>  Pi-  Pn-  then  the  ij  element  of  the  matrix  Q,  after  some  manipulation,  reduces  to 


■ 

V?  f  \  \  fP(0 

k[T'*x** . sm;S,0;I]  =  1  J  ~~  dt 

k=i  i  dQj  de  V(t) 


One  notes  that  each  element  of  represents  the  cross-sensitivity  of  measurement  with  respect  to 
the  response  xsk  of  node  sk. 


The  optimal  sensor  locations  are  then  obtained  by  picking  m  locations  sk,  k  =  1 ,  2 . m,  out 
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of  a  possible  N,  so  that  a  suitable  norm  of  the  matrix  Q  is  maximized  (e,  g.,  the  trace  norm,  etc...), 
[13,  16].  This  may  be  specified  by  the  condition 

max  IIQ[T;si,S2 . sm;S,6;I]  II  .  (90) 

»*£  (1.N) 

3.5.3.  Linear  and  Nonlinear  Systems 

It  is  important  to  note  that  while  nothing  was  said  about  the  nature  of  the  governing 
differential  equations  describing  the  system  models,  the  methodology  presented  up  to  now  is  quite 
valid  for  systems  with  non-linear  differential  equations.  It  should  also  be  noted  that  the  criterion 
developed  by  equation  (89)  does  not  hinge  upon  the  linearity  of  the  system.  The  only  equations 
involved  are  the  measurement  equation  (76)  and  the  relation  (84).  Therefore,  the  methodology 
introduced  herein  may  be  applied  to  the  systems  governed  by  non-linear  differential  equations. 

3.6.  Choice  of  Matrix  Norms 

As  mentioned  above,  Q  is  a  matrix  and  therefore  it  is  necessary  to  use  a  suitable  scalar  norm 
of  it,  IIQII,  to  obtain  an  idea  of  the  information  content,  about  a  parameter  vector  8,  available  from 
sensors  at  one  or  more  locations,  given  the  input,  l(t). 

Various  norms  may  be  used  as  scalar  measures  of  performance.  Some  commonly  used 
norms  are[19]: 

(1)  D  -  optimality:  Minimize  the  determinant  of  Q'1  or  eqivalently  maximize  the  determinant 
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(2)  A  -  Optimality:  Minimize  the  trace  of  Q'1 

(3)  T  -  Optimality:  Maxinize  the  trace  of  Q. 

An  important  advantage  of  D  -  optimality  is  its  invariance  under  scale  changes  in  the  parameters  and 
linear  transformations  of  the  output.  However  T  -  optimality  has  the  advantage  that  the  trace 
operator  is  linear  and  therefore  Trace(Q)  can  be  expressed  as 


When  m  is  large  this  relationship  allows  the  optimal  sensor  locations,  got  by  maximizing 
Trace(Q)  (as  given  by  (90),  to  be  obtained  in  a  simple  sequential  manner.  The  algorithm  to  be  used 
can  be  described  in  the  following  three  steps. 

Step  (1):  For  each  s.  ,  k  =  1,2, .  . ,  N,  determine  q  ; 

* 


Should  r  sensors  be  already  fixed  in  place  at  locations  sk,  k  =  1,  2,..r,  the  best  locations  for 
an  additional  m  sensors  can  be  found  by  including  a  further  step  in  the  above  alogrithm.  After 
performing  Step  (2)  perform  the  following  two  steps. 


Step  (2'):  From  the  sorted  array  obtained  in  Step  (2)  above,  delete  q  ,  k=l,2,..,r  ; 


Step  (3):  The  q  ,  k=l,2, .  . m  locations  corresponding  to  the  largest  values  in  the 
sk 

remaining  sorted  array  yield  the  optimal  sensor  locations  for  the  next  m  sensors. 

Due  to  computational  ease  and  efficiency  with  which  the  Trace  criterion  can  be  used,  and  the 
simplicity  with  which  the  maximization  defined  in  (90)  can  be  carried  out,  in  this  sequel  we  shall 
exclusively  use  the  Trace  criterion.  Comparison  of  the  results  between  A,  T  and  D  -  optimality  will 
be  left  for  a  future  study. 


3.6.1.  Analytical  Interpretation  of  the  Trace  Norm 


In  this  section  we  present  a  more  formal  interpretation  of  the  Trace  norm.  Let  us  introduce  an 
error  criterion 


then 


j  =  E0Y[f(  e,  e  >], 


(93) 


J  =  E 


Yie 


f(e,6)  +iL(§-e)  +  I(0-e)Ti-L($-e 
ae  2  a£ 


(94) 


where  0,  0  and  Y  are  the  true  value  of  the  estimated  parameters,  the  estimate,  and  the  measurement 

A  A 

yielding  the  estimates,  respectively.  Since  f(0,  0  )  is  a  function  of  error  between  the  0  and  0  , 
then  f(0,  0)  =  0.  Hence, 


J  -  E 


(e 

0f  ~ 

—  (9-0) 

+  E 

i(§-6)T  — (§-6)1) 

l  Y  1  8 

La© 

Y  10 

2  a§2  Jj 

(95) 


•  df  df 

However,  if  0  is  close  to  0,  then  — —  =  — .  Using  this  approximation  one  can  write 

00  00 


,“E6(fEv-(S'e)  +  EvjT(§-9,TS<§-e,l) 


(98) 


Notice  that  if  9  is  an  efficient  unbiased  estimator,  then  E(9  )  =  0.  Hence, 


Ev  i  '  ®)  ~  Ev  i  o(9  )  * 


Consequently  one  can  further  simplify  J.  Therefore, 


1  A  T  irf  „ 

J  =  E  E  -(6-6)  — 7(9-6) 
o  y le  2  502 


Noting  that 


Cov$=EYie[($-0)  (0  -  ^  )T  |  0], 


after  some  matrix  manipulation  one  can  write 


I 

J  =  Eft  -Trace  (—Cov§) 

9  2  a§2 


(100) 


To  minimize  the  error  between  the  estimate  §  and  0,  one  would  want  to  minimize  the  right  hand 
side  of  the  above  equation.  If  f  is  quadratic  in  0  then,  the  second  derivative  of  f  with  respect  to  0  in 
(100)  is  a  constant  matrix. 


3.7.  Application  to  a  two-degree-of-freedom  linear  system 


Consider  the  problem  of  finding  the  optimal  sensor  location  in  a  two  mass  system,  modelled 
by  the  two-degree-of-freedom  of  Figure  4,  so  as  to  "best"  identify:  (1)  the  mass  ratio,  A,  of  the 
first  to  second  mass;  (2)  the  stiffness  ratio,  B,  of  the  lower  to  the  upper  spring;  and  (3)  the 
damping  coefficient. 

The  governing  differential  equation  of  motion  can  be  expressed  as 


MX  +  CX  +  KX  =  Wf(t) ;  X(0)  =  0  ,  X(0)  =  0  ,  (101) 

where  X  =  <Xj  x2>T,  C  =  aK,  W  =  <Am  m>T  and  f(t)  is  the  base  excitation.  The  matrices  M  and 
K  are 


(102) 


Let  Sj  denotes  the  lower  mass  location  and  s2  the  upper  mass  location.  The  selection 
between  the  locations  can  be  equated  to  determining  the  one  non-zero  element  of  selection  matrix,  S 
[1x2]’  measurement  H(t)  defined  by  H(t)  =  SX  +  V(t),  where  V(t)  is  assumed  to  be 

stationary  Gaussian  White  Noise  (SGWN)  with  \(/(t)  =  \y0.  As  Xj  denotes  the  lower  mass 
response,  if  S  =  [1  0}  the  lower  mass  is  selected  for  measurement;  if  S  =  [0  1]  the  upper  mass  is 
selected.  The  location  Sj  would  then  be  preferred  over  the  location  s2  for  identifying  the  parameter 


A  (for  instance),  if  QtT.Sj]  >  Q[T,  $2],  where  T  is  the  time  that  the  measurement  is  taken, 


QjCD  =  QjCT,  sj)  =  \ 
Vo 


(103) 


and. 


A  1  (■  dx.  dx, 

Q2(D^Q(T,s2)=-t]  — 

V0oVdA  1 


(104) 


Since  there  is  only  one  parameter  being  estimated  (say  A),  the  Fisher  Information  matrices  reduce 
to  scalar  quantities. 


The  dependence  of  the  OSL  on  various  type  of  the  base  excitations  has  been  studied  by 
considering  various  base  accelerations.  In  all  the  examples  that  follow  it  will  be  asumed  that  the 
parameter  values  corresponding  to  the  MJC,  and  C  matrices  (of  equation  (69))  are  provided  in 
appropriate  and  consistent  units. 

3.7.1.  Example  of  the  2-degree-of-freedom-system  (2-DOFS)  subjected  to 
impulsive  base  excitation 

Consider  the  OSL  problem  for  the  "best"  (minimum  covariance)  identification  of  the 
parameter  A  (given  the  parameter  B  and  a)  of  the  two-mass  system  governed  by  the  differential 
equation  (28)  using  an  impulsive  base  excitation  (f(t)  =  8  (t) ).  Figure  5-A  shows  the  plots  of  the 
ratio  of  the  information  matrices  (T)/Q2(T),  for  T  =  50  seconds,  for  various  values  of  the 
parameters  A  (which  is  to  be  identified)  and,  a*  =  ato0,  where  co0  =  Vk/m  .  Points  on  the  graph 
with  ordinates  greater  than  unity  indicates  the  optimal  location  to  be  the  lower  mass  level  and 
vice-versa.  The  closed  form  solution  for  the  undamped  case  is  provided  in  the  Appendix  A. 
Figure  5-A  indicates  that  the  optimal  location  in  most  cases,  for  the  range  of  A  considered,  is  the 
upper  mass  level.  However,  we  observe  that  for  some  small  values  of  A  and  a*  the  OSL  is  the 
lower  level.  We  note,  interestingly  enough,  that  the  optimal  sensor  location  for  identification  of  A 
actually  depends  not  only  on  the  value  of  B  and  a*  which  are  presumably  known,  but  also  on  the 
value  of  the  parameter  A  which  itself  is  to  be  identified!  Thus  to  be  able  to  ascertain  the 
optimal  sensor  location  some  apriori  assessment  of  A  is  necessary.  This  brings  us  back  to  our 
discussion  of  section  3.2  where  we  emphasized  that  the  optimal  sensor  location  problem  is  looked 


ze< 


Variation  of  the  damped  natural  feguencles  (normalized 


at  herein  in  the  "local"  context. 

Consider  next  the  OSL  for  the  "best"  identification  of  parameter  B  (given  A  and  a).  The 
exact  solution  of  this  problem  is  also  provided  in  Appendix  A.  Figure  5-B  shows  that  the  optimal 
location  for  identification  of  parameter  B,  using  an  impulsive  base  input,  is  again  the  upper  mass 
level  for  the  range  of  B,  however,  and  a*  >  0.05  (  a*  =  aco0,  Uq  =  Vk/m  )  the  trend  appears  to  be 
more  and  more  in  favor  of  the  upper  mass.  This  seems  intuitively  correct  for  as  B  becomes 
larger,  the  lower  mass  becomes  immobile  and  the  OSL  would  be  the  upper  mass. 

3.7.2.  Example  of  the  Two  degrees-of-freedom-system  (2-DOFS)  subjected  to 
sinusoidal  base  excitation 

Consider  the  OSL  for  "best"  identification  of  parameter  A,  B,  and  a  of  the  two-mass  structure 
governed  by  the  differential  equation  (101)  when  f(t)  =  A  Sin(cot). 

Figure  6-A  and  6-B  show  for  convenience  how  the  system  damped  natural  frequencies, 
normalized  with  respect  to  to0  ( co0  =  Vk/m  ) ,  change  with  changes  in  A  and  B.  In  Figure  6-A,  the 
effect  of  A  values  in  the  vicinity  of  A  <  1  on  the  natural  frequencies  is  highly  noticeable.  In  Figure 
6-B,  however,  the  frequency  curves  become  less  sensitive  to  changes  in  B  for  large  enough  values 
of  B. 

Figure  7- A  is  associated  with  the  solution  of  the  OSLP  for  estimating  the  parameter  A  or  B 
using  sinusoidal  base  excitation.  This  figure  shows  that  as  the  normal!  :ed  driving  frequency  y  (y  = 


Variation  of  Q1/Q2  with  y-gj/wq  f°r  different 


co/co0)  varies,  the  OSL  changes.  The  exact  forms  of  the  Fisher  Information  matrices  are  expressed 
in  Appendix  B.  We  observe  that  for  the  case  when  a*  =  0, 


(105) 


when  estimating  the  parameter  A  (or  B).  Figure  7-A  indicates  that  the  OSL  is  always  the  upper 
mass  level  for  y  <  42  and  the  lower  mass  for  y  >  V2  .  The  dependence  of  the  Fisher  value  ratio 
a*  is  also  illustrated.  It  is  of  some  significance  to  also  note  that  the  information  at  both  floors 
vanish  (for  estimating  A)  when  a*  =  0  as  y  — >  1.  Thus  the  region  in  frequency  space  around  y  =  1 
is  not  a  good  region  for  exciting  the  base  of  the  structure  sinusoidally  if  the  aim  is  to  obtain  records 
for  the  identification  of  A.  The  information  from  records  available  at  either  level  is  however  of 
equal  value  (so  far  as  estimating  A  is  concerned)  for  y  =  0  and  y  =  ^2  . 

For  the  estimation  of  B  (given  A  and  a*  =  0)  the  dimensionless  driving  frequency 

Y  =  yj  1  +  —  (derived  in  Appendix  B  )  yields  no  information  on  B  from  records  at  either  of  the 

two  mass  levels.  On  the  other  hand,  the  responses  at  the  two  mass  levels  yields  identical  amounts 
of  information  on  B  at  y  =  0  and  y  =  V2  (for  A  *  1),  as  indicated  by  the  value  of  Q,/Q2  =  1  at  these 
frequencies. 

The  value  of  the  Q1/Q2  =  0  at  y  =  1  is  indicative  of  the  fact  that  the  upper  mass  level  is  a  far 
better  location  for  a  sensor  (when  estimating  B)  and  a*  =  0. 
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White  Noise.  Integration  was  done  over 
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Figure  7-B  indicates  the  OSL  solution  for  the  situation  wherein  a*  (a*  =  aco0)  is  to  be 
estimated.  The  figure  indicates  the  dependence  of  the  OSL  on  the  driving  frequency  y  and  the 
normalized  damping  value,  a*,  for  the  system.  For  the  particular  system  chosen,  for  example, 
when  y  <  0. 16,  the  OSL  is  the  upper  mass  level.  For  values  of  y  >0. 1 6,  the  OSL  depends  on  the 
value  of  the  critical  damping  =  a  con  /2,  where  con  is  the  natural  frequency  of  the  system  (n  =  1 , 
2).  As  the  damping  increases,  the  OSL  moves  towards  the  upper  mass. 


3.7.3.  Example  of  the  2-degree-of-freedom-system  (2-DOFS)  subjected  to  zero 
mean  Gaussian  white  noise  (ZMGWN)  base  excitation 

Consider  the  same  two  mass  structure  governed  by  differential  equation  (101).  Let  f(t)  be 
taken  as  a  sample  of  Gaussian  White  Noise  (GWN).  Figures  8-A,  8-B,  and  8-C  are  associated 
with  the  OSLP  for  estimating.  A,  B,  and  a*  (a*  =  a£00)  respectively.  For  each  case  five  different 
realizations  of  a  GWN  process  were  used.  The  expected  values  of  Qj/Q2  are  indicated  in  each  plot 
with  a  solid  line.  The  dashed  lines  indicate  the  one-standard  deviation,  (1-a),  band  in  which  the 
ratio  of  the  Fisher's  lie.  The  results  indicate  that  for  damped  structural  systems  (a*  >  0), 
the  OSL  in  most  cases  is  the  upper  mass  level.  Hence,  one  would  want  to  locate  the  sensor  on  the 
top  mass  to  obtain  records  yielding  estimation  of  A,  B  and  a*. 

3.8.  Optimal  Sensor  Locations  for  N-Degree-of-Freedom  Linear  Systems 


In  this  sequel  the  solution  to  the  OSLP  will  be  generated  for  multi-degree-of-freedom  linear 
systems.  Consider  the  N-degree-of-freedom,  classically  damped,  dynamic  system  whose 


governing  differential  equation  of  motion  can  be  written  as  : 


MX  +  CX  +  KX  =  F(t),  X(to)  =  X0,  X(to)  =  X0, 


(106) 


where  Xq  and  X0  are  the  given  initial  conditions  for  the  system.  Using  normal  mode  method  the 
response  vector  X(t)  can  be  determined.  Introducing 


X(t)  =  Oq(t) 


(107) 


where  O  is  the  (NxN)  weighted  model  matrix  (transformation  matrix)  and  q(t)  may  be  refered  to 
as  the  N- vector  of  generalized  coordinates  (response  coordinates).  Then, 


ti+2^nconti+Ati  =  <DTF(t),  q(tg)  =  <|>TMX0,  q^)  =  <|>TMX0, 


(108) 


where  the  (N  x  N)  diagonal  matrix  X  (sometimes  known  as  the  generalized  stiffness  matrix)  is 


[  A  ]  =  <DTK<D  =  Diag[  c^,  cOj, .  .  ,  co^,] 


(109) 


The  solution  of  equation  (107)  is  given  as 


TljO)  =  q0  Ujft-t^+fiQ  v.(t  - t^  +  f  h.  (t-x)p  (x)dx 

i  i  * 

*0 


(110) 


where  q0i  and  qa  are  the  initial  conditions  and, 


74 


Ui(t)  =  EXP^oxt) 


Coscod  t+  — -  Sinco^t 

»  co.  » 

a. 

i 


v.(t)  =  - EXP(-^icoit)Sincod  t. 


COj 


hj(t)  =  v,(t), 


and 


Pi(t)  =  <DTF(t),  i  =  1,2,...,N. 


(111) 

(112) 

(113) 

(114) 

(115) 


If  the  OSLP  is  to  be  solved  for  estimation  of  the  parameter  vector  0,  where  0  was 
previously  defined  in  eq.  (70),  one  should  differentiate  equation  (105)  with  respect  to  0.  This 
yields: 

MX0  +  CX0  +  KX0  =  F0(t)  -  (M0X+C0X+K0X);  X@(t  ^  =  0,  X0(t  $  =  0  (116) 


where 
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[M  X  ]  =  {  M  X  I  M  'X  I.  .1  Mo  X) 


(118) 
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Introducing 


yields 
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for  i  =  and 

j  =  1 . L. 

X0  =  <Dz(t) 


z  +  2^0^35  +  Az  =  G(t);  z(tQ)  =  0 ;  zCt^  =  0 


where 


G(t)  =  OT[F0  -  (M0X  +  C0X  +  K0X)]  . 
Equation  (125)  can  further  be  simplified  using  (123)  to 
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(126) 


G(t)  =  <P  [Fq-  (M0Oq  +  C0O>ri  +  K0tDr])  ] 


where  T|  and  T)  can  be  obtained  by  proper  differentiation  of  eq.  (109).  This  may  be  shown  as 
follow 


Tli(t)  =  T|0  w.(t-t0)  +  ii0  Yj  (t  - 10)  +  |  h.  (t  -  x)pi(T)dx 


where 


(£.o).)2 

W.(t)  =-  EXP(-£.(0.t)  CD,  +  1  -  Sincodt, 


(128) 


Yj(t)  =  EXP(-^co.t)  Cosco^ t-l  l  Sincod t  , 


(129) 


hj(t)  =  Y;(t),  and 


(130) 


Pi(t)  =  <DTF(t)  ,  i  =  1,2,...,N  . 


(131) 


i 

ityO  =  Tlo.  W;(t  -  to)  +  ti0  Y.(t  -  t(j)  +  J  h,(t  -  x)Pi(x)dx 


(132) 


where 
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[C0d  +  (4-COj)2]  CoSCC,  t  , 


I  (^°V 

Yj(t)  =  EXP(-tcat)  <  -I-L-  -  cod  Sin©  t 

“i  1  1 


(134) 
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h:  (t)  =  Y.(t) 


p,(t)  =  <DTF(t),  i  =  1,2,...,N. 


(135) 


(136) 


Therefore,  substituting  equations  (127)  and  (132)  into  equation  (126)  gives  G(t).  Consequently 
the  solution  of  equation  (124)  can  be  written  as: 


Zj/O  =J  h.(t  -x)G.j(t)d(x) 


(137) 


where  h;(t)  is  the  same  as  in  eq.  (1 13).  Notice  that  the  initial  conditions  in  eq.  (124)  are  zero.  This 
is  due  to  the  fact  that  the  initial  conditions  of  (106)  are  known  constants. 


If  in  the  eq.  (106)  we  assume  that  [C]  is  expressed  as  a  linear  combination  of  [K]  and  [M], 
then  eq.  (126)  can  further  be  simplified.  Namely, 


C  =  2cxK+2pM 


(138) 


where  a,  and  p  are  known  constants.  Hence  in  equation  (108),  the  percentage  of  damping,  £N, 


can  be  expressed  as: 


=  ctco.  +  -8-  ,  i  =  1,2,...,N 


(139) 


Then  eq.  (126)  may  be  written  as: 


G(t)  =  C>  {F  -[M  <D(ri'+2pTi)  +  K  0(q+2ari)  +  2(p  M+a  K)Oif}  (140) 

0  0  0  0  0 


To  futher  simplify  eq.  (140),  let  us  consider  the  following  three  cases: 


(1)  If  vector  0  contains  @M  only  (estimating  masses),  then 


G(t)  =  O‘[F0-Me<D(Ti+2pft)]  . 


(2)  If,  however,  vector  0  of  eq.  (70)  only  contains  the  subvector  ©K  (estimating  stiffness), 


(141) 


G(t)  =  O  [F0  -  K0<D(Tl+2PTi)]  .  (142) 


(3)  Finally  if  vector  ©  is  reduced  to  subvector  ©c  ,  ©c  =  [  a  (5]T ,  then 


G(t)  =  «DTF  -2Ari  ■  ©TF  -  2In> 

a  .  p 


where  I  is  the  identity  matrix.  It  should  be  realized  that  if  the  input  F(t)  is  independent  of  ©,  then 
F0  would  be  cancelled  all  through  this  discussion. 


Once  the  solution  of  eq.  (116)  is  obtained  using  (123),  the  Fisher  matrices  may  be  evaluated 
using  eq.  (85).  Hence, 
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3.8.1.  Some  useful  results 


Let  us  assume  that  for  a  given  N-degree-of-freedom  dynamic  system  one  seeks  the  solution 
of  the  OSLP  for  the  purpose  of  all  the  system  model  parameters  {0}L.  Let  R=  I  and  ¥(0=  'Fq. 


With  the  equation  governing  the  motion  expressed  (69),  following  the  solution  steps  presented  in 
the  previous  section,  one  can  write: 


'•W1*  -  ^ 


NxL 


(146) 


To  obtain  the  Fisher  Information  Matrix,  we  shall  utilize  equations  (85)  and  (88).  Hence  using  eq. 
(88)  one  can  write: 


0 

0 


I 

*K 


«—  Sj^  -  th  row,  (147) 


0 


and 


H  =  <DTI  . 

*k  \ 

Therefore,  equation  (144),  under  our  assumptions  with  t0=  0,  can  be  written  as: 


(148) 
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where  m  is  the  number  of  sensors  to  be  used.  The  above  equation  can  further  be  expanded  to: 
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k=i  g 
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zH  z. 
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2lh,  *l 


dt  (150a) 


Equation  (150)  is  the  Fisher  Information  matrix  for  the  given  N-degrec-of-freedom  dynamic 
system.  If  0,  is  not  to  be  estimated,  then  the  i-th  row  and  the  i-th  column  of  the  matrix  of  equation 
(150)  would  be  absent.  Therefore,  if  0j  and  03  of  ©  were  to  be  estimated,  then  the  first  and  the 
third  rows  and  columns  in  the  matrix  equation  (150)  would  only  be  presented.  Namely, 
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(150b) 
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Figure  9,  Multi-degree  of  Freedom  Systems  for  Numerical  Studies 


This  result  is  useful  for  computational  purposes. 

3.9.  Applications  to  Multi-Degree-of  Freedom  Systems 

Figure  9  shows  two  different  multi-degree-of  freedom  systems  that  will  be  used  to  illustrate  the 
OSLP  methodology  developed.  The  numerical  results  obtained  will  indicate  the  nature  of  the 
solutions  for  the  OSLP  and  oftentimes  their  non-intuitive  character.  To  illustrate  the  dependence  of 
the  OSL  on  the  nature  of  the  location  (s)  and  types  of  inputs,  two  types  of  excitations  have  been 
used  --  transient  and  impulsive. 

3.9.1.  Fixed-Fixed  System  of  Figure  9(A) 

Figure  9(A)  shows  a  four  degree  of  freedom  system.  The  system  parameters  are:  m,=m2=2; 
m3=m4=l;  kj  =1^=1 00;  k3=75;  k4=50;  k5=50.  The  damping  is  taken  to  be  of  Raleigh  form  (i.e.  C= 
2aM  +  2|3K)  with  a=0.001  and  (3=  0.04.  The  measurement  noise  'P(t)  is  taken  to  be  \|/0. 

Further,  it  is  assumed  that  we  have  the  ability  to  apply  an  impulsive  force  f(t)  (to  any  one  of 
the  masses  n^,  i  =  1, 2,3,4)  whose  impulse 

oo 

1  =  J  f(t)dt  =  10  <151> 

o 

The  parameter  values  provided  are  assumed  to  be  in  appropriate  and  consistent  units.  We  shall 
investigate  the  following  : 
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1)  If  the  impulsive  force  described  above  is  applied  to  one  of  the  masses,  say  mass  m ,  j  e 
(1,4),  then  where  should  we  locate  a  sensor  to  best  identify  one  of  the  stiffnesses  ki ,  i  e  ( 1 ,5)? 

2)  Were  we  required  to  place  more  than  one  sensor  to  identify  k;  how  would  we  find  the 
optimal  locations  ?  Could  we  rank  order  the  locations  1,2,3,4  shown  in  Figure  9(A)  indicating  the 
order  in  which  they  should  be  populated  by  sensors  so  as  to  best  identify  kt  ? 

3)  Can  we  get  an  idea  regarding  the  information  gained  (or  reduced)  by  placing  a  sensor  at 
location  r  as  opposed  to  location  k  (rjc  e  ( 1 ,4))  ? 

4)  Given  that  we  want  to  identify  lq  using  an  impulsive  force  which  can  be  applied  at  one  of 
the  masses,  at  which  mass  should  it  be  applied  and  at  which  locations  should  the  corresponding 
responses  be  measured  for  best  identification? 

5)  What  are  the  answers  to  questions  (1)  through  (4)  above  if  we  want  to  identify  not  just 
one  stiffness  kt  but  a  group  of  them,  say  kj  and  kj  ? 

Figure  10(A)  shows  how  the  information  on  the  stiffeness  parameter  kj  changes  with  time 
for  records  obtained  at  various  locations  when  an  impulsive  force  is  applied  at  mass  trij  with  1=  10 
units.  As  seen  from  the  figure,  the  information  obtained  from  Location  2  (see  Figure  9(A))  is  the 
maximum  and  this  says  that  the  sensor  (for  identification  of  kj),  if  only  one  such  sensor  be 
available,  should  be  placed  at  Location  2.  We  note  that  our  intuitive  idea  of  using  a  sensor  at 
Location  1  would  have  provided  about  36%  less  information  than  the  optimal  sensor  location 
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Figure  10.  Optimal  Sensor  Locations  for 
Impulsive  Force  Applied  to  Location  No.  1 
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obtained.  The  graph  also  gives  the  rank  ordering  for  optimal  sensor  locations  as  {  Location  2, 
Location  3,  Location  1,  Location  4},  Location  2  being  the  best.  Location  4  being  the  worst. 


Figures  10(B),  10(C),  10(D)  and  10(E)  indicate  the  optimal  sensor  locations(OSL)  for 
identification  of  the  parameters  k^,  k3,  k4,  k5  respectively  using  an  impulsive  force  applied  at 
Location  1.  We  note  that  the  OSL  depend  on  the  parameter  that  is  required  to  be  identified.  Also  of 
interest  is  the  fact  that  the  Fisher  Information,  Q(t),  at  each  location  is  a  function  of  time.  Thus 
Figure  10(B)  shows  that  were  to  be  identified  using  simply  a  two  second  length  of  record 
beginning  at  zero  time,  then  Location  1  would  be  the  optimal  location.  However,  the  use  of  a 
longer  duration  of  record  for  identification  of  would  yield  Location  2  as  the  optimal  location  as 
seen  in  Figure  10(B). 

Figures  1 1,  12  and  13  show  similar  results  for  identification  of  the  stiffness  parameters  using 
an  impulsive  force  (with  1  =  10)  applied  at  masses  m^  m3  and  m4  respectively.  We  see  that  the 
extent  of  information  obtained  about  a  parameter  for  the  purposes  of  identifying  it  and  therefore,  in 
general,  the  optimal  sensor  locations,  depend  on  the  location  where  the  force  is  applied.  Thus 
figures  10(A)  and  11(A)  show  that  the  information  about  the  parameter  k,  from  measurements 
taken  at  Location  2  is  about  1.7  times  greater  if  the  impulse  is  applied  at  mass  nij  rather  than  at 
mass  nij.  In  fact  figures  10(A), 11(A),  12(A)  and  (13a)  show  that  to  identify  kj  using  one  sensor, 
the  best  location  for  both  applying  the  impulsive  force  and  for  obtaining  a  measurement  record,  is 
Location  2. 

Figures  14(A)  indicates  results  for  the  situation  where  both  kt  and  k5  are  to  be  simultaneously 
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Impulse  at  Mass  4 


identified  using  noisy  measurements  at  one  or  more  locations  with  an  impulsive  force  (  1=10) 
applied  at  mass  m3.  While  the  identification  of  alone  would  show  the  OSL  to  be  at  Location  2 
(Figure  12(A)),  and  that  of  kj  alone  to  be  Location  5  (Figure  12(E)),  the  OSL  for  simultaneous 
identification  of  both  these  parameters  is  Location  3.  The  Locations  can  be  rank  ordered  as  ( 
Location  3,  Location  2,  Location  4,  and  Location  1),  Location  3  being  the  best.  Should  more 
sensors  be  available,  they  would  then  successively  populate  the  mass  Locations  as  per  this  ordering 
so  that  identification  of  these  two  parameters  can  be  best  carried  out  Figure  14(B)  shows  a  similar 
result  except  that  the  impulse  (I  =10)  is  applied  now  at  mass  m4.  We  observe  that  the  rank  ordering 
of  locations,  as  per  our  Trace  criterion,  has  now  significantly  changed  to  {Location  4,  Location  3, 
Location  2,  Location  1 } ,  Location  4  being  the  best.  Having  answered  the  five  questions  that  were 
posed  above,  we  next  go  on  to  verify  some  of  our  results. 


3.9.2.  Verification  of  Optimal  Sensor  Locations  as  Determined  by  the 
Methodology 

Consider  the  results  depicted  in  Figure  13(D)  where  the  OSL  is  obtained  for  the  situation 
where  an  impulsive  force  (1=10)  is  applied  to  mass  m4  and  identification  of  k4  is  intended.  The 
system  parameters  are  those  chosen  in  Section  3.9.1.  The  correct  value  for  k4is  50  units.  The 
results  show  (Figure  13(D))  that  Location  4  is  far  superior  to  Location  3.  Parameter  identification 
using  the  reccursive  prediction  error  method  (RPEM)  was  carried  out  using  records  obtained  from 
Locations  3  and  4  in  response  to  the  impulsive  force  at  m4.  For  comparison  purposes  the  same 
identification  scheme  (and  computer  program)  was  used  for  records  obtained  from  both  locations. 
The  identification  scheme  was  started  off  with  a  close-by  initial  guess,  namely  k4=40.  The  results 
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of  this  identification  are  shown  in  figures  15(A)  through  15(D)  for  four  different  levels  of 
noise-to-signal  ratios  (indicated  in  the  figure  by  N/S).  To  provide  a  feel  for  the  extent  of  noise 
prevalent  in  the  records  we  show  in  Figure  16  and  Figure  17  the  noisy  displacement  records  used 
for  identification  for  two  noise  levels:  N/S=  0.2  and  N/S  =  0.5.  We  note  that  while  Locations  3  and 
4  provide  about  the  same  accuracy  of  identification  for  N/S  values  less  than  0.1,  as  the  N/S  ratio 
increases,  Location  4  as  predicted  by  Figure  13(D)  is  indeed  superior.  In  fact  identification  of  k4 
can  be  carried  out  with  reasonable  accuracy  even  when  N/S=1.0  as  seen  in  Figure  15(D)  with 
measurements  from  Location  4.  On  the  other  hand  measurements  from  Location  3  cause  the  same 
identification  scheme  to  diverge  for  N/S=0.5  and  N/S=1.0. 

3.9.3.  Fixed-Free  System  of  Figure  9(B) 

We  present  now  an  example  of  OSL  for  a  transient  excitation  provided  at  the  base  of  the 
system  depicted  in  Figure  9(B).  The  system  parameters  are:  m1=m2=2;  m3=m4=l;  m5  =0.5; 
k1=k2=100;  k3=75;  k4=50;  k5=50.  The  damping  is  again  taken  to  be  of  Raleigh  form  (i.e.  C=  2aM 
+  2(3K)  with  a=0.001  and  p=  0.04.  Figure  18  shows  the  base  acceleration  and  the  results  for 
optimal  sensor  location  for  identification  of  parameter  kj.  It  is  interesting  to  note  that  though 
considerations  of  uniqueness  in  identification  would  dictate  Location  1  to  be  optimal, 
considerations  of  identifying  the  parameter  kj  starting  from  a  close-by  estimate  shows  that  Location 
4  is  optimal.  This  example  therefore  brings  out  the  difference  between  "global"  convergence  and 
"local"  convergence  as  discussed  in  Section  3.2. 
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3.10.  The  Kinetic  Energy  Criterion 


Some  investigators  have  proposed  that,  heuristically  speaking,  displacement  sensors  be  located  at 
locations  where  the  kinetic  energy  of  the  system  is  a  maximum.  While  this  may  be  a  somewhat 
intuitive  approach  to  the  problem  we  have  found  that  this  Kinetic  Energy  Criterion  does  not  yield, 
in  general,  the  optimal  sensor  locations.  Firstly,  such  a  criterion  is  not  dependent  on  the  parameters 
that  are  required  to  be  identified  as  any  such  criterion  should.  Secondly,  and  perhaps  more 
importantly,  our  concern  in  locating  sensors  for  best  identification  of  parameters  hinges  around  the 
sensitivity  of  measurements  to  the  parameters  to  be  identified  and  not  on  the  kinetic  energy  of  the 
system.  For  the  sytem  considered  in  Figure  9(A)  the  results  shown  in  Figure  12(E)  indicate  that  the 
optimal  sensor  location  is  at  Location  4.  Figure  19  shows  the  response  of  the  system  at  Locations  3 
and  4  to  an  impulsive  force  (1=10)  applied  at  mass  m3.  Figure  20(B)  shows  that  the  Kinetic 
energy(KE)  at  Location  4  is  lower  than  that  at  Location  3.  Yet  Location  4  is  obtained  as  the  OSL 
from  Figure  12(E).  This  is  explained  by  figure  20(A)  which  shows  that  though  the  KE  at  Location 
4  is  lower  than  at  Loaction  3,  the  sensitivity  (actually,  its  absolute  value)  of  measurement  to 
parameter  k5  is  higher  at  Location  4  than  at  Location  3.  The  results  of  the  KE  criterion  for  the  cases 
depicted  in  figures  14(A)  and  14(B)  are  shown  in  figures  21(A)  and  21(B)  respectively.  Again  we 
see  that  the  KE  criterion  would  lead  to  a  different  and  erroneous  rank  ordering  of  locations  for  the 
simultaneous  identification  of  k,  and  k5. 

3.11.  Conclusions 

We  have  in  this  section  developed  a  methodology  for  optimally  locating  sensors  for 
parameter  identification  using  noisy  measurement  data.  The  methodology  is  predicated  on  starting 


any  such  identification  process  with  a  near-by,  close  initial  guess  estimate  of  the  parameters  to  be 
identified. 

The  optimal  sensor  locations(OSL)  are  shown  to  depend  on:  1)  the  nature  of  the  system  ( the 
structure  of  the  differential  equations),  2)  the  specific  parameter  values  of  the  different  parameters 
in  the  system  model,  3)  the  number  of  sensors  to  be  used,  4)  the  duration  of  time  over  which  the 
identification  is  to  be  carried  out,  5)  the  specific  parameters  to  be  identified,  and  6)  the  nature  and 
location(s)  of  the  input  time  functions  (applied  forces).  A  simple  algorithm,  which  is 
computationally  efficient,  has  been  developed  for  obtaining  the  OSL. 

The  methodology  has  been  applied  to  a  specific  set  of  multi-degree-of-freedom  systems  and 
numerical  results  have  been  obtained.  The  methodology  provides  a  rank  ordering  of  the  locations 
from  best  to  worst.  It  also  answers,  in  a  rational  manner,  where  to  locate  additional  instruments, 
given  that  several  are  already  in  place  in  a  dynamic  system. 

The  methodology  has  been  validated  (though  in  a  limited  sense  only)  by  actually  using  data 
from  various  locations  and  showing  that  those  locations  that  are  predicted  to  be  optimal  by  the 
methodology  do  indeed  provide  tha  best  identification  of  the  parameters  concerned  from  noisy 
measurement  data. 

The  results  have  shown  that  the  heuristically  obtained  Kinetic  Energy  Criterion,  which  is 
sometimes  alluded  to  in  the  literature,  has  little  to  do  with  optimally  locating  sensors  and.  to  that 
extent,  is  inappropriate  for  use  in  developing  methodologies  relevant  to  such  problems. 
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SECTION  4.  PRELIMINARY  RESULTS  IN  DATA  COMPRESSION 


This  section  deals  with  finding  the  optimal  measurement  locations  for  a  structural 
system  modelled  by  a  single-degree-of-freedom  oscillator,  so  that  any  one  of  the 
parameters  to  be  identified  can  be  estimated  with  a  minimum  variance.  The  measurements 
are  assumed  to  be  taken  in  a  noisy  environmment,  and  the  section  addresses  both  linear 
and  nonlinear,  nonhysteretic  systems.  Besides  the  analytical  relations  deduced  for  the 
optimal  measurement  locations,  it  is  found  that,  in  general,  there  may  exist  measurement 
locations  at  which  no  additional  information  on  the  parameter  under  consideration  is 
generated.  For  the  linear  case,  the  optimal  measurement  locations  are  found  to  be 
independent  of  the  system  response  and  the  actual  values  of  the  parameters  to  be  identified. 
They  solely  depend  on  the  nature  of  the  excitation  used  in  the  identification  prodecure. 
Analytical  results  relating  to  the  optimal  measurement  locations  for  minimizing  the  sum  of 
the  variances  of  the  estimates  of  some  of  the  parameters  are  also  provided. 

4.1.  Introduction 

The  identification  of  parameters  in  dynamic  models  for  structural  systems  is  a  field  that 
is  rapidly  gaining  importance.  In  this  section,  we  attempt  to  study  the  optimal  spacing  of 
measurements  for  a  structural  system  modelled  by  a  single-degree-of-freedom  oscillator  so 
that  the  variance  of  one  of  the  parameters  being  identified  is  minimized. 

We  start  with  a  linear  oscillator  and,  using  Fourier  transforms,  derive  a  set  a  linear 
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algebraic  equations.  The  condition  on  the  measurement  frequencies  so  that  the  estimated 
variance  (from  noisy  measurement  data)  of  either  the  mass  parameter,  the  stiffness 
parameter  or  the  damping  parameter  is  minimal  is  derived  analytically.  It  is  also  found  that 
there  may  exist  a  set  of  frequencies  to,  at  which  no  additional  information  on  that  parameter 
is  available,  yielding  no  reduction  in  its  estimated  variance. 

The  determination  of  the  optimal  measurement  frequencies  depends  solely  on  the  nature 
of  the  forcing  function  used  in  the  identification  procedure  and  is  invariant  with  respect  to 
the  actual  values  of  the  parameters  being  estimated. 

A  numerical  example  is  indicated  to  illustrate  the  analytically  obtained  results.  The 
method  is  then  extended  to  find  the  optimal  measurement  times  for  structural  systems 
modelled  by  general  nonlinear,  second-order  differential  equations  which  represent 
memoryless  systems.  It  is  shown,  again,  that  a  set  of  time  points  may  exist  at  which 
measurements,  if  made,  will  not  yield  additional  information  on  the  parameter  of  specific 
concern. 

Though  the  results  obtained  at  this  time  are  pruely  analytical,  it  is  anticipated  that  they 
will  help  in  the  design  of  experiments,  especially  where  data  handling  and  reduction  are  a 
major  cost  concern. 

4.2.  Problem  statement 

Consider  a  structure  modelled  by  a  single-degree-of-freedom  system  oscillator 


subjected  to  an  excitation  force  q(t).  If  x  is  the  displacement  of  such  an  oscillator,  then  its 
equation  of  motion  is 

mx  +  cx  +  kx  =  q(t),  x(0)  =  0,  x(0)  =  0,  (152) 

where  the  parameters  m,  c  and  k  denote  the  mass,  damping  and  stiffness  respectively  and 
are  assumed  to  be  real  numbers.  Taking  fourier  transforms,  this  yields 

-mo^Xfco)  +  iccoX(co)  +  kX(co)  =  Q(co),  (153) 

so  that 

l/X(co)  =  (-mco2  +  ico>  +  k)/Q(to),  to  e  In,  (1 54) 

where  we  shall  assume  that  the  division  on  both  sides  of  equation  (153)  is  possible  i.e., 
there  exists  an  open  interval,  In,  such  that  for  all  to  e  In,  X  (co)  and  Q  (co)  are  not 
identically  zero. 

We  shall  assume  that  the  parameters  m,  k  and  c  need  to  be  identified,  and  we  shall 
direct  our  interest  to  finding  if  there  exists  a  set  of  frequencies  co;  i  =  1,  2,...  such  that  the 
variance  of  any  one  of  the  desired  parameters  (i.e.,  m,  or  k  or  c)  can  be  minimized  by 
using  the  data  i.e.,  X  (co)  and  Q  (co)  at  those  specific  frequencies. 


Relation  (154)  can  be  rewritten,  after  separating  the  real  and  imaginary  parts,  as 


U(co) 

-co  f(co) 

-cog(co) 

f(CO) 

UJ 

V(co)_ 

— 

-co2g(co) 

cof(co) 

g(co)_ 

c 

k 

+ 

cpv(t») 

(155) 


where  we  have  represented 

l/X(co)  =  U(co)  +  iV(co), 

1/Q(g>)=  f(co)  +  ig(co),  (156) 

and  the  measurements  U  (co)  and  V  (to)  are  corrupted  by  measurement  noise  Ey  (co)  and 
£v(co).  Defining 


[a/co)] 


A 

*  a  (co) 


A 


I-co^co)  -cog(co)  f(co)]T 


(157a) 


and 


[b/co)]  =  b(co)  -  [-co2g(co)  cof(co)  g(co)]T  (157b) 

for  each  value  of  co  =  co;,  i  =  1,  2,...,  N,  equation  (155)  can  be  written  and  the  BLUE 
estimator  obtained.  This  can  be  expressed  by  the  relation 


z  =  H9  +  e, 


(158) 


z  =  [U(co1)V((o1)U(co2)V(co2)-U(0)n)V(con)]T, 
H  =  [a(cOj)  b^)  a(co2)  b(o>2)-a(con)  b(can)]T, 
e  =  [e^)  eytco^  ev(co2)-eu(o)n)  ev(con)]T. 


(159) 


Assume  that  the  error  vector  e  has  the  statistic 


E[e]=0  and  R  =  E(eeT)  =  diag (a-i2.<J]2,a22 . ar^.a^  ),  (160) 

where  a{  =  a  (co^,  i  =  1,  2,...,  n.  The  covariance  of  the  BLUE  estimate  of  the  vector  0 
becomes  [16] 

P  =  (HtR'1H)*1.  (161) 

where  P  is  a  3  x  3  matrix  whose  diagonal  elements  Pu,  P22  and  P33  are  the  variances  in 
the  estimate  of  m,  c  and  k  respectively. 

4.3.  Optimal  choice  of  frequencies,  co*. 

We  shall  now  attempt  to  choose  the  frequencies  a\e  I  n  in  such  a  way  that  the  i  th 
element  of  P,  P-,  is  minimized.  To  that  end  we  first  differentiate  relation  (161)  with 
respect  to  to^  to  yield 


2  -t  2 


dP 


=  -Pi 


^R'H  +  HViL-HViLR'H 


9to, 


3co,. 


do). 


P, 


(162) 
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yf 

V, 


A 


I 


^V* 


i 


v-v.v>» 


W  =-2p^r‘hp  +  phV 

IT  K  J  . 


Assuming  that  a(to),  b(u)  and  a(to)  are  continuous  functions  in  I 


OH7^)  =  [0...0  afo^)  bfc^)  0...0], 


so  that 


—  {p[a(tok)aT(cok)  +  b(cuk)bT(cok)]p} 

^ \  o[ 

+  {p[a(o)k)aT(tok)  +  b(cok)bT(cok)]  p} 


On  expanding,  equation  (164a)  becomes 


•  Z  {p,a/“k,a.(Mk)p» + } 
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0(01) 

4- - — 

o(o)*) 


(164b) 


The  condition  that  Pn  be  extremal  then  yields 


Xv/v 

m' 

j* 

✓ — \ 

£ 

_ t 

+ 

Xwv 

X  wv 

j 

L  *  J 

i 

S 

0(0)*) 


o(o>*)  [ 


XvJ  ♦  Xw 


v  j 


\  j 


=  0. 


(165) 


LEMMA  1.  If  a  real  to*  e  Iq  exists  such  that  f(to*)  and  g(o)*)  are  not  zero,  and  for  any  r 
e  [1,2,3], 


2.W-o. 


(166) 


then. 


2,  Psrbs(o)*)  =  0,  for  that  value  of  r, 


(167) 


and  vice  versa. 


Proof.  Let  us  say  that  we  have  a  sequence  of  CD's,  CDj,  co2,...,  co^  where  for 

relation  (166)  is  valid.  Using  relations  (159)  and  (161)  we  have. 


3_!  _  ^  a(cos)aT(CDs)  +  b(cos)bT(cos) 


so  that 


N  4  2.  , 

V  tt5a  (^s) 

s=i  a2 


N  2  2., 

y  msa 


s=i  a 


N  2  2., 

y 

s=i  a2 


N  2  2., 

£  ((°s> 


y  aVop 


s=i  a 


where 


a2(cos)  =  f^co.)  +  g2(co ). 


The  determinant  A  of  P-1  then  becomes 


(168) 


(169) 


vr»-p»T»T». 


I 


>  ">  --  J- 


A-I^VCs). 


N  2/  \  N  4  2..  22.. 

Y  a2(£0s)  Y'  C0sa  (C0s)  ^  “s®  (°V 

2  Ami  2  2 


^  2  2 
s=i  a  s=i  a 


,  (170) 


which  by  the  Cauchy-Schwartz  inequality  is  always  >  0,  as  it  should  be,  since  P  is  a 


covariance  matrix.  Thus  the  matrix  P  now  becomes 


AB  0  B 


P=  0  CA-B1  0  /A, 

2 

n**  A  A  n 


0  BC 


where 


V1  a2(a>s)  'Y' 

A  =  X-T-’  B=  X 


r  ©  a2(o)s) 


5—  1  O 


s=i  a 


N  4  2.  . 

c  =  ^  {05a  ^ 


s=i  a 


A  =  B(AC  -  B2). 


Noting  relation  (157),  condition  (166)  becomes 


f(o>k)[-G)k  A  +  B]  =  0,  r=l 

-g(tO|t)  [to^CA  -  B2)]  =  0,  r  =  2, 

f(cok)[-  B  +  C]  =  0,  r  =  3, 


(171a) 


(171b) 


(171c) 


(172) 


(173a) 


Relation  (167)  yields 


g(<*\)  [-  «>*  A  +  B]  =  0,  r  =  1, 
fCcc^)  [co*  (CA  -  B2)]  =  0,  r  =  2,  (173b) 

g(co^)  [- ©£  B  +  C]  =  0,  r=3. 

If  g(c*\)  and  fta^)  *  0  then  two  sets  of  equations  become  identical,  and  the  result 
follows.  ■ 


LEMMA  2. 


and 


£Pnbl  . 

IP0a,  f<“k>  ’ 


,  Pljaj  *  0,  r  *  1,3 


ZPA  _  fK>  t 

IPrjaj  «<«V)  ’ 


*  0,  r  =  2. 


(174) 


(175) 


Proof.  The  proof  follows  from  equations  (173a)  and  (173b).  ■ 

THEOREM  1.  For  a  given  forcing  function  Q(co)  and  any  r  €  (1,2,3),  there  may  exist 
frequencies  ci^  such  that  the  inclusion  of  data  at  those  frequencies  does  not  yield  any 
improvement  in  the  variance  of  our  estimate  of  parameter  r.  Specifically,  when  coke  IQ 
satisfies  equation  (166),  to  =  ci^  is  such  a  frequency. 


V  '.f’J'TXW  •TO'ICW  V  v.v  v.v.v.  v.  --.  --.  w. 


.  *  -  *  . 


Proof.  Let  us  imagine  that  the  measurements  at  the  frequencies  tOj.  co2,...,cok 4  have 
been  made  and  that  with  each  measurement,  the  covariance  matrix  P  is  updated.  After 
making  the  k  th  measurement  at  o  =  ci^,  the  updated  covariance  matrix  becomes  [16] 

P*  =  P  -  P  Hj  [Rk  +  HkP  Hj  ]  **  HkP,  (176) 


where. 


P  denotes  the  covariance  before  the  measurement  at  o^, 

P  +  denotes  the  covariance  after  the  measurement  at  co^ 

Hk  denotes  [a,b]T  evaluated  at  o^,  and 
Rk  =  a|  diag  (1,1). 

Relation  (176)  can  be  rewritten  using  the  notation 

L  ■  P  Hj  (177) 

as 

P*  =  P  -  L  [Rk  +  HjjP  H^L1.  (178) 

Using  relations  (157)  and  (159),  we  have 
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L  = 


1 V, 

IP„b, 

1  Vi 

IP2ibi 

iP3jbi 

(179) 


If  relation  (166)  is  valid  for  some  r,  then 

X  p„ a,  -  X  p„  b«v  <180> 

and  by  relation  (179) 

Lrj  =  0’  j  =  1,2  for  that  r.  (181) 

Consequently,  from  equation  (178)  we  find 

p;,=P;(  j  =  1,2,3.  (182) 

We  note  in  passing,  using  relations  (180)  and  (182), 

IP+a.  =  ZP+b=0.  (183) 

n  j  n  j 

COROLLARY  1.  For  values  of  tok  for  which  X  Prjaj(o)k)  =  0  (for  r  =  1,3),  the 
variance  of  the  estimate  of  the  r  th  variable  as  well  as  its  covariance  are  unaffected  by  the 
measurement  co  =  co^ 
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Proof.  The  result  follows  directly  from  relation  (182).  ■ 

We  note  that  the  frequencies  oij.,  which  do  not  contain  any  further  information  (for 
any  given  r),  do  not  depend  on  the  parameter  values  m,  c  and  k.  They  are  only  governed 
by  the  nature  of  the  forcing  function  Q(o>)  and  can  be  calculated  before  even  the 
measurements  are  made.  They  also  do  not  depend  on  the  measured  responses. 

COROLLARY  2a.  If  g(cok)  *  0  and  f(o>k)  *  0,  there  exists  no  (Dke  In  except 
possibly  =  0  for  which 

I  P2jaJ(cak)  =  0’  084) 

provided  the  covariance  matrix  is  nonsingular. 

Proof.  Referring  to  equations  (173)  for  r  -2  and  noting  that  AC  -  B2  *  0  because  the 
determinant  of  the  covariance  matrix  is  nonzero,  the  equations  can  only  be  satisfied  by 
=  0. 

COROLLARY  2b.  If  g(o\)  *  0  and  f(o\)  *  0,  there  always  exists  an,  to^  e 
(0,  <~)  ,  for  which  LP^fa^)  =  0,  r  =  1,3. 

Proof.  We  shall  show  the  case  for  r  =  1.  Using  relation  (173a),  coj[  =  B/A.  From 
equation  (172),  B  >  0  and  A  >  0.  Thus  =  (B/A)1/2. 

The  proof  for  r  =  3  follows  along  similar  lines.  ■ 


THEOREM  2.  The  optimal  locations  for  the  measurements  0  =  0^,  o^e  I  n,  if  they 
exist  at  all,  which  minimize  the  variance  in  the  estimates  Pfft  satisfy  the  following  relations: 


g(cO 

1 P  a  (ox ) - -  IP  b  (col)  - 

0  J  %  f(a^)  ^ 


■ 

\2 

O(03k) 

1  + 

g(Wk)| 

<J(C\) 

I  P  a  =  0,  r  =  1,3  (185a) 

rj  J 


when  f^)  *  0.  When  f(03^)  =  0,  they  satisfy  the  relation 

o(£a )  ^ 

P rb  (a^)  -  — -  V P„b(03k)  =  0,  r  =  1 .3,  (185b) 

*  o(03k)  , 

Proof.  Using  equation  (165)  and  Lemma  2,  the  result  follows.  A  similar  result  can  be 
written  for  the  cases  when  g(o>k)  *  0  and  g(0)k)  =  0  respectively. 

Equations  (185a)  and  (185b)  express  the  criteria  for  finding  observation  points 
such  that  the  mass  m  or  the  stiffness  k  can  be  optimally  identified.  We  note  that  the 
optimal  location  of  the  kth  observation  point  depends  in  general  upon  the  location  of  all 
the  previous  observation  points  as  contained  in  Prj  and  Psr . 


THEOREM  3.  The  optimal  locations  for  the  measurements  03  =  03^  03ke  I  fl,  if  they 
exist  at  all,  which  minimize  the  variance  of  the  damping  parameter  c,  satisfy  the  relation 


(186) 


(0,0(0^) 

<j((Dk) 


lg2(e\)  +  f2^)] 


^  —  [g2(C0)  +  f^co)]  I 
2  do 


Proof.  Noting  relation  (165)  and  Lemma  2,  the  result  follows.  We  assume  that  the 
covariance  matrix  is  strictly  positive  definite. 

Theorem  3  states  that  to  optimally  locate  the  k  th  measurement,  relation  (186)  needs  to 
be  satisfied.  We  note  that  in  this  case  the  optimal  location  of  the  k  th  measurement  does 
not  depend  on  the  locations  of  the  proceeding  measurements  and  is  purely  controlled  by  the 
nature  of  the  graphs  of  f(oj)  and  g(co). 

Next  let  us  consider  the  problem  of  minimizing  the  sum  of  the  variances  of  X  out  of 
the  3  parameters.  Let  s(,  i  =  1,...,  X  be  these  X  parameters.  Let  A  be  the  zero  matrix 

whose  nth  diagonal  element  is  unity  if  n  *  s{,  i  *  1 . X,  Sj  e  (1,  2,  3).  We  then  have  the 

followng  result. 

THEOREM  4.  The  frequencies  In  which  make 

2 

I  P„.  ^  S3,  (187) 

«=i  *J 


extremal  are  given  by  the  relation 


a  (o>k)PA 


Pa(o>, ) - Pa(co. ) 

+  b  (cok)PA 

Pb(C0  ) - - 

a(cok) 

o(cok) 

=  0, 


where  A  is  the  selection  matrix  as  defined  above. 


Proof.  The  extremal  condition  is  given  by 


I 

i  =  1 


dP 


SS 

■ii  =0. 


dco* 


Using  relation  (164a)  this  becomes 


from  which  the  results  follows.  ■ 


COROLLARY  3.  The  extremal  values  of  trace(P)  are  given  by  the  relation 


(188) 


(189) 


090) 


aT(wk)  P2a(cok)  +  bT(cok)  P2b(cok)  =  0 


•  091) 


when  a  =  0.  This  corresponds  to  R  =  aQ  Diag(  1,1,  .  .  ,1)  . 


Proof.  For  this  case  A  =  I  and  the  result  follows.  ■ 


COROLLARY  4.  If  either  or  gCci^)  is  nonzero,  then  there  exists  no  >  0, 

for  which  the  vectors  P  a((Dj.)  and  P  fcKco^)  equal  zero,  provided  P  is  nonsingular. 


Proof.  Noting  that  the  vector  Pa  is  proportional  to 


-  ca^A]  ^ 

'  glto^-AC] 

^f(o)k)[C  -  0)2B]  J 


(192) 


where  A,  B  and  C  are  defined  in  relation  (171-172),  the  result  follows.  The  result  for  P  b 
is  along  the  same  lines. 

COROLLARY  5.  If  X  >  1,  there  exists  no  frequency  o^,  o^e  IQ,  such  that  for  flo^) 
and  g(wk)  nonzero  and  P  nonsingular 


AP  a  =  0. 


(193) 


Proof.  The  result  follows  from  Corollary  4  and  relations  (173). 


We  note  therefore  that,  in  general,  the  data  at  each  frequency  provide  information  on  m 


and/or  k. 


4.4.  Numerical  example. 


JO,  t  <  0, 

ta,(I)  =  L*,  t  >  0,  p  >  0 


Then,  f(co)  =  (3  and  g(co)  =  co.  The  vectors  a  and  b  become 


(194) 


a  =  [-co2P  -co2  P]T, 


b  =  [-co3  cop  co]1 


(195) 


Let  us  assume  that  measurements  at  co  =  C0j,  t  have  been  made  and  the  next 

measurement  is  to  be  taken  at  cok.  Let  P'  denote  the  covariance  matrix  at  the  end  of  the  first 
(k  -  1)  measurements,  and  C0j  €  [1,  °°).  Let  o(co;)  =  a0,  Vi,  so  that  c(co)  h  0. 

If  cOj,.  is  such  that  X  Pj^co^  =  0,  then  no  improvement  in  the  variance  of  the  stiffness 
estimate  can  be  expected  by  obtaining  the  additional  measurement  at  co  =  cok  This 
condition  for  the  forcing  function  (194),  after  some  algebra,  implies 


cok  =  [C,/B']1/2, 


(196) 


iCvOTO  .a.  k.*  O*  v  O  O  «  ■  •  * 


c=  +cos)  20(1  B'  =  +0i? 

S=1  S=1 


(197) 


For  illustration,  assume  that  k  =  4  and  [3=1.  If  the  first  three  measurements  are 
taken  at  0)  =  1  rad/sec.,  £0  =  2  rad/sec.,  co  =  3  rad/  sec.,  and  o  =  a0,  j  =  1,2, 3, 4,  then 


.00648  0  .04274 

P '  =  o\  0  .0089  0 

.04274  0  .3405 


B’  =  1 12,  C  =  892,  and  relation  (196)  gives  co4  =  2.82  rad/sec.  Thus  a  measurement  at  to 
=  2.82  rad/sec.  will  yield  no  improvement  of  the  covariance  of  the  stiffness.  In  fact,  from 
Corollary  1  we  know  that  P|(=  P3  j,  j  =  1,2,3.  Again  if  o^is  such  that  XP^a^o^)  =  0, 
then  the  new  measurement  will  provide  no  additional  information  on  the  mass  parameter. 
After  some  algebra,  we  find  that  this  relation  gives 


(^  =  (B 7  A’)1/2, 


(199) 


where  A’  =  z£=j(p2  +  ct^2),  and  B’  is  as  defined  before.  We  then  have  P*j  =  Pj.  , 

j  =  1,2,3.  For  the  example  taken,  A’  =  17,  and  the  value  of  satisfying  (199)  is  2.56 
rad/sec. 


For  optimally  locating  the  measurement  0)  =  cok  so  that  Pu  is  minimized,  equation 


ssv*. 


(185a)  needs  to  be  satisfied.  This  yields 


-2(32  ^  a2(cos)  +  3g)2  ^ a2(co  )  +  ^ co2a («)  =  0, 


(200) 


which  simplifies  to 


2co*  +  co2 (4(3 2  +  3A')  +  2A'{32  -  B3  +  2[32  =  0. 


(201) 


For  cd4s  [1,°°)  in  our  example,  relation  (194)  gives  a  value  of  1.148  rads/sec. 


4.5.  Extensions  to  some  nonlinear  sdof-systems. 


In  this  section  we  extend  the  results  obtained  dealing  with  the  parameter  identification  of 


a  linear  system  to  nonlinear  systems  that  can  be  described  by  one  degree  of  freedom.  The 


proofs  of  all  the  results  follow  suit  from  those  of  the  previous  section  and  therefore  have 


been  omitted.  Here  we  shall  work  directly  in  the  time  domain. 


Consider  a  structure  modelled  by  a  nonlinear  differential  equation 


nrx  +  f(x,  -x)  =  q(t),  x(0)  =  -x(0)  =  0, 


(202) 


where  within  a  certain  range'of  response  the  nonlinear  term  can  be  approximated  by 


V*  yv^  %  J*yCTTgi?'V^rVVVTr  it*  v  ’/  "J*  '/•  ’l 

•///.*  A  •  -\  A.  /.  /.  -  ./  -  %  N 


f(x,x)=  XPMxn(t)i”(t). 


(203) 


The  system  is  assumed  to  be  nonhysteretic.  Let  us  say  that  the  response  of  the  system 
is  measured  at  times  t  =  tk ,  k  =  1,2,...,  t  e  (0,  T),  and  that  the  aim  is  to  locate  the  instants 
t  =  tk  when  measurements  should  be  made  so  that  the  data  collected  thereat  would  yield  the 
minimal  variance  of  any  one  of  the  parameters  m  or  (3^  whose  accurate  identificaiton  is 
required.  We  shall  assume  that  x,  x  and  x  and  their  derivatives  are  continuous  functions 
of  time  for  t  e  (0,  T). 


If  the  complete  time  histories  x,  x  and  x  were  available  (actually  they  are  not),  then 
equation  (202)  could  be  rewritten  as 


mzj(t)  +  CjZjO)  +  CjZjO)  +  —  +  cLz  L  +1(0  =  q(0. 


(204) 


where  the  z's  correspond  to  the  corresponding  time  functions.  L  =  m  x  n  and  the 
coefficients  Pmn  are  assembled  into  a  one-dimensional  array  c.  Identification  of  the 
parameter  vector  0  =  [m  cT]T  under  noisy  measurement  conditions  would  lead  to  the 
relations 


q  =  H9  +  e, 


(205) 


where 


Vs'1.  ’<  V’  ’V* 


(206) 


q  =  [q(tj)  q(t2)  -  q(tN)]x, 

H  =  [z(t1)z(t2)...z(tL+1)]T,  z(t)  =  [Z|(t)  Z2  (t) ...  zL+1(t)]T, 
e  =  the  zero-mean  white  measurement  noise. 


The  covariance  of  the  estimate  can  be  written,  as  before,  as  the  (1  +  mn)  x  (1  +  mn)  matrix 

P  =  (HTR^H)-1,  (207) 


where  R  is  the  noise  covariance  matrix  and  is  taken  to  be  diag  (of  -  cf  ,...,(*£ ).  Once  again 
the  extremal  condition  for  P  with  respect  to  a  measurement  instant  t^  can  be  expressed,  as 
before,  by  the  relation 


s  - 

Iw 

j 

S 

(208) 


where  z-  is  the  j  th  element  of  z.  (We  note  that  this  result  follows  directly  from  equation 
(165)  if  we  set  b  to  0,  a  to  z,  and  co^  to  tk.)  We  then  have  the  following  result. 

Theorem  5.  For  a  given  forcing  function  and  any  r  6  {1,2,...,  mn  +  1 },  there  may  exist 
times  tj.  such  that  the  inclusion  of  data  at  those  times  does  not  yield  any  improvement  in  the 
variance  P^  of  our  estimate  of  the  r  th  parameter.  Specifically,  when  for  some  tk  e  (0,T) 
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mn=l 

L  W° 

j=l 


(209) 


is  satisfied,  then  t  =  ^  is  such  an  instant. 

Proof.  The  proof  follows  along  the  lines  of  Theorem  1.  ■ 

Corollary  6.  For  values  of  tk  for  which  lPrjzj  (^  )  =  0,  the  variance  of  the  estimate 
of  the  r  th  variable  as  well  as  its  cross-covariance  is  unaffected  by  the  measruement  at  t 

-*k- 

Proof.  As  before,  the  condition  implies  that 

pq  =  Pq  .  j  =  l^,...mn  +  1.  (210) 

The  proof  follows,  as  did  that  of  Corollary  1,  from  the  proof  of  Theorem  5.  ■ 


Theorem  6.  The  optimal  measurement  times  tk,  tk  €  (0,  T),  if  they  exist  all,  which 
minimize  the  variance  satisfy  the  following  relation 


Proof.  The  proof  follows  along  the  lines  of  Theorem  2.1 


THEOREM  7.  The  measurement  times  ^  e  (0,  T)  which  cause 


P,  s  ,  l  £  mn  +  1, 


to  be  extremal  are  given  by  the  relation 


T  g(l  )  j 

z'O^PAPz^)  — i-zVgPAPz^) 
<y(tk) 


when  A  is  the  zero  matrix  whose  n  th  diagonal  element  is  1  if  n  =  sif  i  =  1,...,  Jt,  s 
(1 . mn  +  1). 


Proof.  The  proof  is  exactly  along  the  lines  of  Theorem  4. 


Theorem  8.  If  z  (tk )  is  such  that  for  a  given  set  s;,  i  =  1,  2,...,  1;  X  <,  mn  +  1, 


APz(t.)  =  0 


i 

then  there  exist  time  t,  €  (0,T)  such  that  I  ,  P  remains  unatfeeted.  In  fact.  P 

k  1=1  ss  rs 

1  1  1 

r  =  1,  2,...,  mn  +  1,  i  =  1,2,...,  Z  remain  unchanged  by  the  new  measurement  at  time  tk  . 

Proof.  The  proof  follows  along  the  same  lines  at  that  of  Theorem  1.  ■ 

4.6.  Remarks  and  Conclusions. 

In  this  section  we  have  tried  to  understand  the  optimal  measurement  strategy  for 
identifying  the  parameters  of  a  single-degree-of-ffeedom  dynamic  system. 

For  a  linear  system,  we  have  shown  that,  in  general,  data  acquired  at  all  frequencies,  in 
the  interval  cok  e  Ifl,  do  not  equally  enhance  our  knowledge  of  the  parameters  being 
estimated.  Specifically,  one  can  often,  given  a  data  stream  collected  at  frequencies  00^1  = 
1,2,...,  k  -  1,  forecast  the  next  frequency  at  which  data  collection  would  be  maximally 
beneficial  to  obtaining  a  more  confident  estimate  of  any  one  of  the  desired  parameters. 
Likewise,  one  can  predict  the  frequency  at  which  data  collection  would  have  no  influence 
on  improving  the  uncertainty  in  our  estimate  of  any  desired  parameter.  It  is  shown  that  for 
{ co  €  IQ  /  f(co),  g(co)  *  0 ) ,  data  at  all  the  frequencies  carry  information  about  the  damping 
parameter  c.  Also  as  opposed  to  the  optimal  measurement  location  for  identification  of 
m  and  k,  which  do  depend  on  the  previous  measurement  locations  to,,  co2,--  cok the 
optimal  locations  for  the  identification  of  c  do  not  depend  on  the  locations  of  the 
measurement  stream.  They  are  solely  controlled  by  the  nature  of  the  forcing  functions 
used  in  the  identification  procedure. 


Irrespective  of  which  parameter  is  being  identified,  the  optimal  measurement  locations 
do  not  depend  on  the  values  of  the  parameters.  It  is  noted  that  the  solution  of  equations 
(185)  and  (186)  may  not  exist  for  any  co  belonging  to  the  open  interval  IQ.  For  such 
situations  the  optimal  locations  would  have  to  be  chosen  as  the  end  point(s)  of  the  interval. 
Also,  it  is  observed  that  the  optimal  measurement  locations  do  not  depend  on  the  system 
output  x(t).  Thus  the  optimal  locations  can  be  calculated  a  priori  to  obtaining  the 
measurement  stream.  A  numerical  example  has  been  included  to  illustrate  the  analytical 
results  obtained.  Furthermore,  analytical  results  relating  to  minimization  of  one  or  some 
sum  of  the  variances  of  the  estimates  are  also  provided. 

The  results  are  extended  to  nonhysteretic  nonlinear  oscillators  when  time  data  of 
displacement,  velocity  and  acceleration  can  be  obtained  for  a  given  time -dependent  forcing 
function.  It  is  shown  that  there  may  exist  certain  times  ^  at  which  measurements  provide 
no  additional  information  regarding  any  one  of  the  parameters  desired  to  be  most  accurately 
estimated.  Similarly  there  exist  times  tk  at  which  measurements  provide  maximal 
information  about  a  desired  parameter.  The  relations  that  these  times  tk  satisfy  in  each  of 
the  two  cases  have  been  analytically  deduced. 

The  results  of  this  section,  it  is  hoped,  will  shed  light  on  the  manner  in  which 
experimentation  can  be  performed  so  that  the  amount  of  data  handling  and  reduction 
required  could  perhaps  be  significantly  decreased  in  the  dynamic  testing  of  structural  and 
mechanical  systems. 
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SECTION  5.  TRADEOFFS  BETWEEN  IDENTIFICATION  AND  CONTROL 

IN  DYNAMIC  SYSTEMS 

5.1.  INTRODUCTION 

The  control  of  large  flexible  structures  is  an  area  that  has  attracted  considerable  interest  in 
recent  years  from  both  the  professional  and  the  research  community  17],  [21-24],  This  interest  is 
primarily  motivated  by  the  need  to  precisely  control  flexible  structures  in  various  developing  fields 
of  modem  technology.  In  the  area  of  earthquake  engineering  the  reduction  of  structural  response 
may  be  necessary  to  reduce  internal  stresses  caused  by  dynamic  loads  thereby  reducing  the  damage 
potential  and  increasing  the  useful  life  of  the  structure  [25].  In  space  applications  the  availability  of 
Space  Shuttle  to  transport  large  payloads  into  orbit  at  reasonable  costs  presents  an  opportunity  for 
large  systems  to  perform  new  missions  in  space.  However,  because  of  launch  weight  and  volume 
constraints  these  structures  are  generally  very  flexible  and  pose  new  challenges  in  all  aspects  of 
control  such  as  attitude  control  and  maneuvering,  precision  pointing,  vibration  attenuation  and 
structural  and  shape  control  [7]. 

A  necessary  prelude  to  the  effective  control  of  a  structure,  is  a  knowledge  of  its 

characteristics  and  properties.  In  other  words,  one  needs  to  have  information  about  the  structural 

system  so  that  adequate  control  algorithms  can  be  devised.  This  has  led  to  a  considerable  interest  in 

the  identification  of  structures  subjected  to  dynamic  loads[8-10].  For  structural  systems  that  are 

described  by  parametric  models,  this  involves  knowledge  of  the  nature  of  the  governing  differential 

% 
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equations  and  knowledge  of  the  values  of  the  parameters  that  are  involved,  or  at  least  know  ledge  of 
the  bounds  within  which  the  parameters  lie.  Clearly,  the  better  the  system  is  identified  (the  smaller 
the  bounding  intervals  within  which  the  parameters  are  known  to  lie),  the  more  finely  tuned  the 
controller  can  be  made,  so  that  for  a  given  amount  of  available  control  energy,  the  control  would  be 
more  efficient.  The  less  knowledge  we  have  about  the  structural  system  the  more  robust  the 
controller  needs  to  be,  and,  in  general,  the  less  efficient  the  control.  Thus  heuristically  speaking 
there  exists  a  duality  between  the  concepts  of  identification  and  control,  because  1)  robust 
controllers  may  require  reduced  efforts  at  identification  (for  purposes  of  control),  and,  2)  increased 
efforts  at  identification  may  require  less  robust  and  more  efficient  controllers.  However  the 
tradeoffs  between  identification  and  control,  from  a  practical  standpoint,  are  still  usually  difficult  to 
assess  quantitatively.  Little  work  has  been  reported  to  date  in  this  area  of  quantitative  cost-benefit 
analysis  between  these  two  dual  concepts. 

In  this  section  we  formulate  the  trade-off  problem  between  identification  and  control,  and 
study  in  a  quantitative  manner  their  duality  through  the  use  of  the  intermediary  concept  of  an 
optimal  input.  Thus  the  section  attempts  to  answer  the  following  question:  given  that  the  optimal 
input  time  function  is  to  have  a  certain  prescribed  energy,  how  does  it  change  in  character  and  in  its 
effectiveness  as  one  changes  the  objective  criterion  from  one  that  emphasizes  control  to  one  that 
emphasizes  identification?  While  the  analytical  work  presented  here  has  been  motivated  by  our  need 
to  control  flexible  structural  systems  the  methodology  developed  and  the  results  obtained  are 
applicable  to  all  systems  governed  by  ordinary  differential  equations;  thus  it  applies  to  systems 
commonly  met  with  in  chemical,  civil,  electrical,  and  mechanical  engineering.  Some  simple 
numerical  examples  are  provided  to  indicate  the  quantitative  nature  of  the  results  and  provide  a  feel 
for  them.  The  results  in  these  examples  show  that  significant  trade-offs  exist  between  identification 
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and  control  and  that  for  the  same  amount  of  energy  in  the  input  signal,  the  emphasis  on  control 
could  lead  to  very  high  covariances  of  the  parameter  estimates.  Similarly  inputs  that  are  optimal  for 
identification  could  yield  responses  whose  mean  squared  values  may  be  several  times  those 
obtained  for  inputs  that  yield  optima]  control. 


5.2.  PROBLEM  FORMULATION 


Consider  a  dynamic  system  modelled  by  the  first  order  set  of  differential  equations 


x(t)=  F,  x  +  f 
z(t)  =  H,  x(t)  +  v(t) 


(215) 


(216) 


where  x  is  an  nXl  state  vector,  f  is  an  mXl  control  vector,  z  is  an  rXI  measurement  vector  and  the 
nXl  initial  condition  vector,  x0,  is  given.  We  shall  assume  that  the  measurement  noise  is 
representable  as  a  zero  mean  Gaussian  White  Noise  process  so  that 


E[  v(t)  ]  =  0,  and, 

E(  v(t)  v(x)  ]  =  R,  5  ( t  -  x  ) . 


(217) 

(218) 


Let  the  vector  of  unknown  parameters  in  the  system  modelled  by  equations  (215)  and  (216)  be 
given  by  the  pXl  vector  0.  Let  us  assume  that  the  identification  is  carried  out  with  an  efficient 

a 

unbiased  estimator  so  that  the  covariance  of  the  estimate  of  0  namely,  0  is  provided  by  the  inverse 
of  the  Fisher  Information  Matrix[13].  Hence, 


Cov[  0  ]  =  M"1. 


(219) 


The  matrices  Ft  and  Gj  are  taken  to  be  functions  of,  in  general,  the  parameter  vector  0.  The 

optimal  input  for  identification  of  the  parameter  vector  0  is  then  sought  such  that  a  suitable  norm 
related  to  the  matrix  M  is  maximized  or  minimized.  Different  measures  of  performance  related  to  M 
have  been  used  in  the  literature  [  19]: 

1.  A  -  Optimality,  where  Tr(  M*1)  is  minimized; 

2.  D  -  Optimality,  where  the  determinant  of  M'1  is  minimized;  and, 

3.  E  -  Optimality,  where  the  maximum  eigenvalue  of  M*1  is  minimized. 

In  this  section,  for  expository  purposes,  we  shall  use  the  criterion  for  obtaining  the  optimal  inputs 
for  identification  as  the  maximizadon  of  the  Trace  (W 1/2  M  W1/2)  where  W  is  a  suitable  positive 
definite  weighting  matrix.  Thus  the  criterion  for  obtaining  the  optimal  input  for  parameter 
identification  is  taken  to  be 

T 

J,  =  J  Trace]  h[ R',1  H,  Vp  I  dt  (220) 

o 

where 

Vp  =  Xp  W1/2  ,  (221) 

and  the  matrix  Xp  is  given  by, 

d  x. 

[X  ]  =  — -  (222) 

pi°  ae 

j 

In  addition  to  the  objective  function  generated  by  our  need  for  identification,  the  objective  function 
required  to  be  maximized  for  control  is, 


Here  is  a  symmetric  positive  definite,  nXn  weighting  matrix.  This  then  yields  the  composite 
objective  function  which  is  required  to  be  maximized  as 


J 


T  T 

-(ct/2)  J{xT  Q(x)  dt  +  (p/2)  J*Trace  {  yj  hJ  R,1  H,  yp  }  ,  (224) 

o  o 


where  a  nd  P  are  positive  scalars.  Clearly,  when  a  »  P,  finding  f(t)  to  maximize  J  is  tantamount 
to  finding  the  optimal  control  for  the  system  (2 15)-(2 16),  while  when  P  »  a  ,  the  f(t)  that 
maximizes  J  is  simply  the  optimal  input  for  identification  of  the  pXl  parameter  vector  0.  In 
particular,  when  a  =  0,  and  p  =1,  the  optimal  input  for  'best'  identification  is  obtained;  when  a  = 
1,  and  P  =  0,  the  optimal  input  for  'best'  control  is  obtained.  Denoting  the  nXl  vector 


l 


(225) 


and  assuming  that  the  matrix  W  is  diagonal,  so  that, 

W  =  Diag(Wj,  w2, . Wp)  (226) 


we  can  generate  an  augmented  n(p  +  1)  vector. 
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with 


-  _  £!l 

e>  a  e. 


G  = 
9; 


a  g  j 
Te~ 


,  i  =  1,2,  .  .  .,  p. 


(231) 


We  note  in  passing  that  the  stability  of  the  equation  set  (228)  is  controlled  by  the  stability  of  the 
equation  set  (215)  since  the  eigenvalues  of  the  matrix  F  are  identical  to  those  of  the  matrix  Fj,  as 

seen  in  equation  (229),  except  for  the  increased  multiplicities.  The  objective  function  (224)  can  now 
be  rewritten,  after  some  algebra,  as 

T  T 

J  =  -  (ot/2)  J  yTQ  yT  dt  +  (p/2)  J  yT  HT  R'1  H  y  dt  (232) 

o  o 

where,  the  matrices  Q,  H,  and  R*^  are  the  block  diagonal  matrices  given  by 
Q  =  Diag{  Q|,  O,  0 . 0,0}, 

H  -  Diag{0,  Hj,  Hv  Hx . Hj },  and,  (233) 

r1  =  Diag{  r;1  ,  r;1  ,  r;1 . r;1  } 

% 

Thus  the  objective  function  needs  to  be  maximized  under  the  constraint  equations  (228)  and  the 
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energy  constraint 


T 

J  fTf  dt  =  E,  (234) 

o 

where  the  parameter  E  is  given  &  priori. 


5.3.  Determination  of  Optimal  Inputs  for  Simultaneous  Identification  and  Control 

Having  formulated  the  problem  for  constrained  maximization,  we  next  use  the  standard 
Lagrange  multiplier  method  to  obtain  the  function  f(t)  which  maximizes  the  objective  (232)  subject 

to  (228)  and  (234).  Using  the  Lagrange  multipliers  X(t)  and  v(t)  we  therefore  obtain  the  augmented 
objective  function  to  be 


T  T  T 

J  =  -  J  {yTQy}  dt  +  &  J  {yTHTR1Hy}  dt  +  j\T(t)  {  Fy  +  Gf  -  y}dt  + 
0  0  0 


(235a) 


Here  we  augmented  the  state  vector  by  the  variable 


WO  «  J  fT(t)  f(t)  dt  (235b) 

o 

where  N  =  n(p+l)  and  we  use  the  additional  Lagrange  multiplier  rj(t).  By  doing  this  we  can  satisfy 
the  equality  constraint  (234)  without  having  to  resort  to  the  usual  trial  and  error  procedure. 
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Taking  the  first  variation,  we  obtain  the  following  set  of  equations  for  the  variables  y(t),  ^.(t) 


and  T|(t)  (see  Appendix  C): 


y(t)  —  Fy  =  -  VXft)  ,  y(0)  =  yQ 

71 

X(t)  +  FT  X(t)  =  -  a  {Q}y  +  (3  (HT  R'1  H  }y  ,  X(T)  =  0,  (236) 

yN+1=fTf=ixT(t)vx(t),  yN+1(0)  =  o,  yN+1(T)  =  E 
71 

ii(t)  =  o 

where  V  =  (GGT).  We  note  that  the  equation  set  (236)  constitutes  a  nonlinear  two  point  boundary 
value  problem  containing  2[np  +  n+  1]  first-order  differential  equations.  The  optimal  input  vector, 
fit),  is  obtained  through  the  solution  of  this  two  point  boundary  value  problem  using  the  relation: 

fit)  =-  —  GTX(t)  (237) 

71 

It  is  interesting  to  note  that  had  we  used  the  objective  function 

T 

Jc'  =  [  J  xT  Qj  X  dt  ]  1  (238) 

o 


instead  of  Jc  in  equation  (223)  we  would  obtain, 


J' 


J2. 

2 


(239) 


where  Jj  is  given  in  equation  (220),  and  the  relative  weightings  of  the  contributions  of  the  control 


and  the  identification  objectives  are  denoted  by  a  and  P'.  Following  the  same  procedure  as  before 
and  using  the  augmented  vector  y,  maximization  of  this  objective  function,  J',  along  with  the 
constraints  (228)  and  (234)  would  then  yield  the  following  set  of  equations: 

y(t)  -  Fy(t)  =  i  V&t)  ,  y(0)  =  yQ  , 

T1 

X  (t)  +  FT  X(t)  =  -  <x  (Q)y  +  (p'A)  {  HT  R1  H  }y  ,  X(T)  =  0  ,  (240) 

yN+1  =  —  i  T(t)  V  X(t),  y  (0)  =  o,  yN+1(T)  =  e 
Tl2 


T)  (t)  =  0  , 


where,  A  is  the  positive  quantity  defined  by 

T  T 

A  =  [  JyTQ  ydt]2  =  [  J xTQ,x  dt  ]2  (241) 

o  o 

The  optimal  input  is  obtained  from  the  relation 


f(t)=  -1  GT  X(t).  (242) 

Tl 

Comparing  equations  (236)  and  (240)  we  observe  that  the  only  difference  that  arises  in  the  use  of 

the  objective  function  (239)  instead  of  (232)  is  the  effective  change  in  weighting  parameter  p. 
Setting 
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P’A  =  p  the  equation  set  (240)  becomes  identical  with  the  set  (236).  The  lagrange  multipliers  X(t) 

of  (236)  and  X  of  (240)  are  related  by  X  =  XA;  similarly,  T|  =  r\A.  It  should  be  noted  that  the  two 
equations  become  identical  only  when  the  function  y(t)  in  (241)  corresponds  to  the  response  for  the 

optimal  input  f(t),  i.e.,  the  solution,  y,  of  the  set  (236).  With  this  rescaling  of  the  parameter  p’,  the 
optimal  input  f(t)  is  identical  for  the  two  objective  functions  J  and  J'  of  (224)  and  (238).  Having 
thus  shown  the  quasi-equivalence  of  the  two  objective  functions  (224)  and  (239)  through  this 
rescaling,  in  this  sequel,  we  shall  illustrate  our  results  by  using  the  objective  function  in  the  form  of 
equation  (224)  which  leads  to  the  boundary  problem  described  in  equation  (236). 

This  two-point  boundary  value  problem  cam  be  numerically  solved  in  various  ways.  An 
extensive  literature  on  numerical  techniques  for  solving  such  problems  is  available[26].  Among  the 
methods  most  commonly  used  are  multiple  shooting  techniques[26]  with  Newton-Raphson 
iterations,  and  the  Kalaba  Method[27]  where  the  two  point  boundary  value  problem  is  converted  to 
an  equivalent  Cauchy  initial  value  problem.  In  this  sequel,  the  equation  set  (236)  is  solved  using  the 
multiple  shooting  technique  with  Newton  iteration. 


5.4.  Illustrative  Example 


To  exemplify  the  concepts  developed,  let  us  consider  a  system  modelled  by  a  single 
degree-of-freedom  oscillator  described  by  the  differential  equations 


Xj(0)=  a0 ,  x^O)  =  b0 


(243) 


where  f(t)  is  the  optimal  input  to  be  applied.  Denoting  x  =  [x  j  X2]T,  xk  =  dx/3k,  and  xd  =  dx/dc ,  ] 

1 

the  objective  function  is  taken  to  be  ) 
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J(T)  =  -aJc(T)  +  bJk(T)  +  d  Jd(T)  (244) 


where, 

T 

Jc(T)  =  J  x,  dt  .  (245) 

o 

T 

Jk(T)  =  J  \  dt  .  (246) 

o 

and, 

T 

JdCT)  =  J  xd2  dt  ,  (247) 

o 

The  weighting  factors  a,  b,  and  d  are  taken  to  be  non-negative.  This  may  be  thought  of  as  being 
produced  by  choosing  a  =  a  ,  p  =  1,  W=  Diag  (b  ,  d  ),  Q]  =  Diag(l,  0)  and  H,  =  [1  01  and  the 

scalar  =  a2  =  1  in  our  general  formulation  of  equations  (224)  -  (232).  The  six  component 
vector  y  is  then  given  by 

■ 
x 

**  ,  (248) 

A/2 


and  the  matrix  F  becomes 


Since  the  vector  G  =  [  0  1  0  0  0  0],  the  vector  { WX/  rj)  has  only  one  non-zero  component  namely, 

X.2(t)/"n,  where  X2(t)  is  the  second  element  of  the  vector  X(t  ).  Figure  22  shows  some  of  the 

numerical  results  for  the  following  parameters  values  (which  we  shall  assume  are  taken  in 
consistent  units): 

b=0,  T=  5,  E  =  2.25;  and,  (250) 

k=  50,  c=2,  x, (0)  =  0,  >4(0)  =10.  (251) 

These  parameters  thus  look  at  a  single  degree-of-freedom  oscillator  whose  spring  constant  is  50 
units,  and  whose  viscous  damping  is  2  units.  This  yields  a  system  which  has  an  undamped  natural 
frequency  of  vibration  of  about  7  radians/sec  and  a  percentage  of  critical  damping  of  about  15%.  It 
is  subjected  to  an  initial  velocity  of  10  units.  The  aim  is  to  study  the  trade-off  between  1) 
identifying  (in  equation  (243))  the  damping  parameter,  c,  in  the  best  possible  way,  and  2) 
controlling  the  system  so  that  its  mean  square  response  over  the  time  period  T  is  a  minimum,  given 
that  an  \iput  (forcing  function)  of  5  units  duration  with  an  energy  of  up  to  2.25  units  is  to  be  used. 
The  ,wo  point  boundary  value  problem  posed  in  equation  set  (236)  is  numerically  solved  using  the 
standard  multiple  shooting  technique!  15],  The  local  error  tolerence  during  integration  of  the 
differential  equations  and  the  permitted  error  in  the  satisfaction  of  the  boundary  conditions  are  each 
set  to  10" 4.  The  responses  of  the  system  together  with  the  optimal  inputs  as  obtained  from  equation 
(237)  are  shown  for  the  two  extreme  cases:  1)  a=0,  d=l,  corresponding  to  the  optimal  input 


ENERGY  »  30.0 


(b) 


requires  for  identification  of  the  damping  parameter  c,  and,  2)  a=l,  d=0,  corresponding  to  the 
optimal  input  required  for  minimizing  the  response.  As  seen  from  Figure  22(B),  the  optimal  inputs 
required  for  ’best’  identification  and  for  'best'  control  (the  term  ’best’  is  used  in  in  terms  of  the  cost 
function  (244)  utilized)  are  widely  different  form  each  other.  In  fact  they  are  seen  to  be,  for  the 
entire  duration  over  which  they  last,  almost  exactly  out  of  phase.  Differences  in  the  response  of  the 
system  to  combined  influence  of  the  initial  velocity  and  the  forcing  functions  obtained  for  the  two 
cases  are  shown  in  Figure  22(A).  Figure  22(C)  shows  the  Fisher  Information  Matrices  for 
damping,  which  in  this  case  are  scalars,  namely  Jd(t),  for  the  abovementioned  two  extreme  cases, 
as  a  function  of  time,  t.  The  difference  between  these  at  T=  5  is  about  55%.  Alternatively  put,  the 
input  forcing  function,  which  controls  the  system  response  maximally,  causes  a  response  which  is 
only  about  55%  as  informative  about  the  system  parameter  c  as  that  caused  by  a  forcing  function 
that  is  designed  to  maximally  provide  information  about  the  parameter,  c.  The  manner  in  which  the 
integral  of  the  response  quantity  squared,  Jc(t),  changes  with  time  for  the  two  cases  mentioned 

above  is  shown  in  Figure  22(D).  As  seen,  at  T  =  5,  the  optimal  control  input  is  about  35%  more 
effective  in  reducing  the  mean  square  response  than  the  input  which  optimally  determines  the 
parameter  c.  Figure  22(E)  provides  the  sensitivity  of  the  response  to  the  damping  parameter  (at 
c=2)  as  a  function  of  time. 

Figure  23(A)  shows  the  manner  in  which  the  Fisher  Information  matrices  Jd(T)  and  Jk(T) 

change  for  various  values  of  the  ratio,  y  =  [  {a  Jc(T))/{d  Jd(T))]  when  the  Fisher  values  are 

normalized  to  unity.  It  is  to  be  noted  that  the  optimal  input  when  a=0  corresponds  to  the  that 
required  for  'best'  estimation  of  the  parameter  c  in  equation  (243).  Figure  23(B)  shows  the  effect  of 
changing  the  available  control  energy,  E,  from  a  value  of  2.25  to  30  keeping  all  other  parameters 
the  same.  From  a  loss  of  information  in  the  parameter  c  of  55%  in  the  case  of  E=2.25,  the  loss  in 
information  when  E  =  30  jumps  to  about  450%.  Similarly,  the  extent  to  which  the  system's 
performance  can  be  controlled  deteriorates  by  about  a  factor  of  3  if  one  aims  at  purely  identification 

instead  of  control.  The  parameter  rj  which  is  calculated  for  each  objective  function  ratio,  y,  is 
shown  in  Figure  24.  As  mentioned  in  the  formulation,  this  quantity  is  simultaneously  solved  for,  in 
the  set  (236),  thereby  eliminating  the  need  to  find  its  value  by  trial  and  error.  Had  this  not  been 
done  a  very  high  computational  expense  would  have  been  incurred  to  ensure  that  the  energy 


constraint  is  satisfied.  Noting  that  the  inverse  of  Jd(T),  for  an  efficient  unbiased  estimator,  is  the 
covariance  of  the  estimate  of  the  parameter  c,  Figures  25(A)  and  25(B)  provide  the  trade  off 

between  control  and  identification.  As  seen  in  Figure  25(A),  in  going  from  y=  0  to  y  =  4,  JC(T),  the 
mean  square  response,  falls  off  by  about  35%;  similarly  the  covariance  of  the  estimate  of  c 

increases  by  about  55%  as  y  varies  over  the  same  interval.  For  larger  values  of  the  input  energy, 
Figure  25(B)  shows  that  significant  reductions  in  Jc  and  significant  increases  in  the  covariance  of 
the  parameter  estimates  can  occur. 

5.5.  Conclusions 

In  this  section  we  have  presented  an  approach  to  quantifying  the  tradeoff  between  the  tasks 
of  control  low  development  and  plant  identification.  To  the  best  of  our  knowledge  this  tradeoff  has 
never  been  analyzed  quantitatively  before.  The  problem  is  formulated  in  the  context  of  optimal 
control  and  optimal  identification  through  the  intermediary  concept  of  an  optimal  input.  A  suitable 
objective  function  is  chosen  so  that  the  emphasis  from  control  to  identification  can  be  changed  in  a 
continuous  manner.  It  is  shown  that  the  duality  between  identification  and  control  can  be  quantified 
by  determining  optimal  inputs,  which  have  a  specified  amount  of  energy,  and  which  minimize  the 
objective  function.  Augmenting  the  state  by  an  additional  variable  allows  a  simultaneous  solution 
of  the  optimization  problem  together  with  the  energy  constraint.  Using  variational  calculus  this 
leads  to  a  two  point  boundary  value  problem  that  is  nonlinear  due  to  the  introduction  of  the  energy 
constraint.  The  boundary  value  problem  is  solved  numerically  using  the  multiple  shooting 
technique  and  Newton-Raphson  iterations. 

A  numerical  example,  which  deals  with  control  and  identification  of  the  parameters  of  a 
single  degree-of-freedom  oscillator,  is  used  to  illustrate  the  concepts  involved.  It  is  shown  that 
improved  control  leads  to  serious  deterioration  in  the  covariance  of  the  parameter  estimation  and 
vice-versa.  In  general,  as  the  energy  of  the  input  increases  the  tradeoffs  between  identification  and 
control  are  shown  to  become  more  and  more  intense. 

The  above  example  shows  the  potential  of  the  approach  introduced  in  this  section.  The 
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APPENDIX  A 


If  in  the  equation  (101)  f(t)  =  5  (t),  the  expressions  for  Qj  (T)  and  Q2  (T)  when  estimating  A 
(a  =  0)  may  be  written  as: 

Qi  (T)  =  pf  a1  +  q,2  bj  +  p|  aj  +  q|  b2  +  2p1q1c1  +  2p,p2d  +  p^g 

+p2q1h+2q1q2^+2p2q2c2  (A.l) 

and 

=  P 1 2aj+q  j  2bj  +p  2  2a2+q'  2  2b2+V t  q' x  c  x +V  i  P 2  d 

+p  iV  22g+p  2q  th+2qiq  2  *+2p  2q  2  c2 

(A.2) 

where 
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also 
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Let  us  define 
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Let  us  further  define 
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Substituing  equations  (A.3  -  A.23)  into  equations  (A.l)  and  (A.2)  yields  the  exact  form  of  Q,(T) 
and  Q2(T)  for  the  solution  of  the  OSLP  for  determining  the  parameter  A. 

Similarly  if  one  seeks  the  solution  of  the  OSLP  for  determining  the  parameter  B  (under  the 
same  impulsive  base  input)  for  the  system  governed  by  equation  (101)  and  (a  =  0),  one  may  write, 
after  some  algebraic  manipulations: 
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where  a4,  bit  and  c;  (i  e  [1, 2] )  are  given  by  equations  (A.3  -  A. 6);  and  d,  g,  h  and  i  are  given  by 
equations  (A. 9  -  A.  12).  Also 
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Substituting  equations  (A. 26  -  A. 33)  into  equations  (A. 24)  and  (A. .25)  yields  the  exact  solution  of 
Q,(T)  and  Q2(T)  for  determining  the  parameter  B. 
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APPENDIX  B 


a)  Herein  we  present  the  closed  form  solutions  for  Qj  and  Q2  associated  with  the  identification 
of  (1)  parameter  A,  (2)  parameter  B  and  (3)  parameter  a  from  response  data  obtained  using  a 
sinusoidal  base  motion. 

The  solution  to  equation  (101)  when  f(t)  =  cQSin(wt)  (a  >  0)  can  be  shown  as 
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(1)  If  the  OSL  for  the  "best"  estimation  of  parameter  A  is  sought  (using  the  approximation 


— - ) ;  i  =  [1,2]),  Q.(T)  may  be  written  as, 
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where. 
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The  quantities  ^  and  ^  in  the  above  equations  can  be  expressed  as: 
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for  i  e  [1,  2]. 
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Similarly, 
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(2)  If  the  OSL  for  the  "best"  estimation  of  parameter  B  is  sought,  Q,(T)  may  be  written  as 
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—  Cos  (cot  +  2 \\i )  Sin(cot)] 
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Similarly 
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(b)  An  easier  approach  to  finding  the  OSL  for  the  response  measurement  to  identify  A  or  B  can 
be  used  for  the  case  a  =  0.  Writing  the  response  X  =  YSin(ot)  we  have, 


Y  =  [Q] 
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i -m 
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where 


[«]  = 


(B+l)k  -  Amo2 


-k 


-k 


k  -  mo2 


(B.44) 


Hence  we  have,  after  some  simplification, 
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We  observe  that  though  the  ratio  Qi(T)/Q2(T)  tends  to  zero  as  y  tends  to  unity,  indicating  that  the 
upper  story  is  the  prefered  sensor  location,  the  values  of  Qj(T)  and  Q2(T)  becomes  vanishingly 
small  for  y  =  1 .  Also  Qj(T)  =  Q2(T)  for  y  =  0  and  y  =  V2  . 

Similarly,  the  expressions  for  Qj(T)  and  Q2(T)  associated  with  the  identification  of  B  can  be 
expressed  as: 
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APPENDIX  C 


Taking  the  first  variation  of  equation  (235)  we  obtain 


T  t  T 
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+  J{  XT  +  XtF}  8y  dt  +  J  XtG  8f  dt  +Jri(t)8fTf  dt  -  —  8yN+,  lQ 
o  o  o  2 

T 
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Collecting  terms  in  8y,  8f  and  Sq,  we  get  the  relations: 
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